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PREFACE 

This  is  the  third  in  a  series  of  Technical  Notes  (910-  )  entitled  "Self-Study  Manual 
on  Optical  Radiation  Measurements."   It  contains  the  sixth  chapter  of  that  Manual.   Inad- 
vertently, this  chapter  has  been  completed  before  chapters  4  and  5,  which  are  nearly  ready 
for  publication  in  NBS  TN  910-2.   Chapter  4,  "More  on  the  Distribution  of  Optical  Radiation 
with  Respect  to  Position  and  Direction,"  deals  with  all  of  the  radiometric  quantities  of 
Table  Al-3  in  Appendix  A  of  NBS  TN  910-1  that  were  not  covered  in  detail  in  the  first  three 
chapters.   Chapter  5,  "The  Measurement  Equation,"  introduces  the  central  topic  that  ties 
together  our  entire  treatment  of  optical  radiation  measurements.   The  measurement-equation 
concept  is  introduced  by  deriving  the  appropriate  measurement  equations  for  three  different 
illustrative  problems.   All  of  these  chapters  are  still  in  Part  I — Concepts,  but  our  plans 
call  for  us  next  to  turn  to  one  or  two  chapters  in  Part  III — Applications.   We  feel  that  we 
now  have  established  enough  of  the  fundamentals  to  be  able  to  usefully  deal  with  some 
specific  applications,  the  topic  in  which  we  realize  many  readers  will  be  most  interested. 

For  background  information  on  the  whole  project  and  on  the  plans  for  the  "Self-Study 
Manual"  (SSM) ,  we  reproduce  here  (immediately  following)  the  Preface  to  the  first  Technical 
Note  910-1,  of  March  1976.   Although  the  exact  details  of  the  outline  and  plans  for  the 
individual  chapters  are  still  developing,  our  principal  aims  and  overall  plans  are  still 
as  set  forth  in  that  Preface. 

Here,  in  chapter  6  of  Part  I,  we  are  primarily  concerned  with  establishing  the  con- 
cepts needed  to  deal  with  polarization  in  relation  to  radiometry  (the  measurement  of 
optical  radiation) .   There  already  exist  excellent  treatments  of  polarization  that  are 
concerned  primarily  with  polarization  phenomena,  and  we  have  not  tried  to  compete  with  them. 
This  is  not  a  treatise  on  polarization  phenomena,  as  such,  but  only  on  as  much  about  them 
as  is  needed  to  deal  adequately  with  the  measurement  of  power  propagated  in  the  form  of 
optical  radiation  —  with  radiometry.   To  our  knowledge,  this  is  a  topic  that  has  not 
previously  been  treated  in  the  literature. 

We  are  grateful  to  many  for  valuable  comments  and  criticisms  received,  both  informally 
as  well  as  in  formal  reviews,  concerning  all  phases  of  the  SSM.   In  connection  with  this 
polarization  chapter,  we  are  particularly  indebted  to  Drs.  Jean  M.  and  Harold  E.  Bennett 
of  the  Michelson  Laboratory,  Naval  Weapons  Center,  China  Lake,  California,  and 
Dr.  Elio  Passaglia  of  NBS  for  helpful  discussions  and  critical  reviews.   Mrs.  Betty  Castle 
has  again  done  an  outstanding  job  of  typing  a  difficult  text. 

Fred  E.  Nicodemus,  Editor 

Henry  J.  Kostkowski,  Chief, 
Optical  Radiation  Section 

PREFACE  to  NBS  TN  910-1 

This  is  the  initial  publication  of  a  new  series  of  Technical  Notes  (910)  entitled 
"Self-Study  Manual  on  Optical  Radiation  Measurements."   It  qontains  the  first  three  chap- 
ters of  this  Manual.   Additional  chapters  will  be  published,  similarly,  as  they  are  com- 
pleted.  The  Manual  is  being  written  by  the  Optical  Radiation  Section  of  NBS.   In  addition 
to  writing  some  of  the  chapters,  themselves,  Fred  E.  Nicodemus  is  the  Editor  of  the  Manual 
and  Henry  J.  Kostkowski,  Chief  of  the  Section,  heads  the  overall  project. 

In  recent  years,  the  economic  and  social  impact  of  radiometric  measurements  (including 
photometric  measurements)  has  increased  significantly.   Such  measurements  are  required  in 
the  manufacture  of  cameras,  color  TV's,  copying  machines,  and  solid-state  lamps  (LED's). 
Ultraviolet  radiation  is  being  used  extensively  for  the  polymerization  of  industrial  coat- 
ings, and  regulatory  agencies  are  concerned  with  its  effects  on  the  eyes  and  skin  of 
workers.   On  the  other  hand,  phototherapy  is  usually  the  preferred  method  for  the  treatment 
of  jaundice  in  the  newborn.   Considerable  attention  is  being  given  to  the  widespread  utili- 
zation of  solar  energy.   These  are  just  a  few  examples  of  present  day  applications  of  opti- 
cal radiation.   Most  of  these  applications  would  benefit  from  simple  measurements  of  one  to 
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a  few  per  cent  uncertainty  and,  in  some  cases,  such  accuracies  are  almost  essential.   But 
this  is  rarely  possible.   Measurements  by  different  instruments  or  techniques  commonly  dis- 
agree by   10%   to   50%,   and  resolving  these  discrepancies  is  time-consuming  and  costly. 

There  are  two  major  reasons  for  the  large  discrepancies  that  occur.   One  is  that  opti- 
cal radiation  is  one  of  the  most  difficult  physical  quantities  to  measure  accurately. 
Radiant  power  varies  with  the  radiation  parameters  of  position,  direction,  wavelength,  time, 
and  polarization.   The  responsivity  of  most  radiometers  also  varies  with  these  same  radia- 
tion parameters  and  with  a  number  of  environmental  and  instrumental  parameters,  as  well. 
Thus,  the  accurate  measurement  of  optical  radiation  is  a  difficult  multi-dimensional  prob- 
lem.  The  second  reason  is  that,  in  addition  to  this  inherent  difficulty,  there  are  few 
measurement  experts  available.   Most  of  the  people  wanting  to  make  optical  radiation 
measurements  have  not  been  trained  to  do  so.   Few  schools  have  had  programs  in  this  area 
and  tutorial  and  reference  material  that  can  be  used  for  self-study  is  only  partially  avail- 
able, is  scattered  throughout  the  literature,  and  is  generally  inadequate.   Our  purpose  in 
preparing  this  Self-Study  Manual  is  to  make  that  information  readily  accessible  in  one  place 
and  in  systematic,  understandable  form. 

The  idea  of  producing  such  a  manual  at  NBS  was  developed  by  one  of  us  (HJK)  in  the 
latter  part  of  1973.   Detailed  planning  got  under  way  in  the  summer  of  1974  when  a  full- 
time  editor  (FEN)  was  appointed.   The  two  of  us  worked  together  for  about  one  year  devel- 
oping an  approach  and  format  while  writing  and  rewriting  several  drafts  of  the  first  few 
chapters.   These  are  particularly  important  because  they  will  serve  as  a  model  for  the  rest 
of  the  Manual.   During  this  period,  a  draft  text  for  the  first  four  chapters  was  distrib- 
uted, along  with  a  questionnaire,  for  comment  and  criticism  to  some  200  individuals  repre- 
senting virtually  every  technical  area  interested  in  the  Manual.   About  50  replies  were 
received,  varying  widely  in  the  reactions  and  suggestions  expressed.   Detailed  discussions 
were  also  held  with  key  individuals,  including  most  of  the  Section  staff,  particularly 
those  that  will  be  writing  some  of  the  later  chapters.   In  spite  of  the  very  wide  range  of 
opinions  encountered,  all  of  this  feedback  has  provided  valuable  guidance  for  the  final 
decisions  about  objectives,  content,  style,  level  of  presentation,  etc. 

In  particular,  we  have  been  able  to  arrive  at  a  clear  solution  to  difficult  questions 
about  the  level  of  presentation.   Both  of  us  started  out  with  the  firm  conviction  that, 
with  enough  time  and  effort,  we  should  be  able  to  present  the  subject  so  that  readers  with 
the  equivalent  of  just  elementary  college  mathematics  and  science  could  easily  follow  it. 
That  conviction  was  based  on  our  experience  of  success  in  explaining  the  subtleties  of 
radiometric  measurements  to  technicians  at  that  level.   What  we  failed  to  consider,  how- 
ever, was  that,  in  making  such  explanations  to  individuals  we  always  were  able  to  relate 
what  we  said  to  the  particular  background  and  immediate  problem  of  the  individual.   That's 
just  not  possible  in  a  text  intended  for  broad  use  by  workers  in  astronomy,  mechanical 
heat-transfer  engineering,  illumination  engineering,  photometry,  meteorology,  photo-biology 
and  photo-chemistry,  optical  pyrometry,  remote  sensing,  military  infrared  applications,  etc. 
To  deal  directly  and  explicitly  with  each  individual's  problems  in  a  cook-book  approach 
would  require  an  impossibly  large  and  unwieldy  text.   So  we  must  fall  back  on  general 
principles  which  immediately  and  unavoidably  require  more  knowledge  and  familiarity  with 
science  and  mathematics,  at  the  level  of  a  bachelor's  degree  in  some  branch  of  science  or 
engineering,  or  the  equivalent  in  other  training  and  experience. 

In  its  present  form,  the  Manual  is  a  definitive  tutorial  treatment  of  the  subject  that 
is  complete  enough  for  self  instruction.   This  is  what  is  meant  by  the  phrase  "self-study" 
in  the  title.   The  Manual  does  not  contain  explicitly  programmed  learning  steps  as  that 
phrase  sometimes  denotes.   In  addition,  through  detailed,  yet  concise,  chapter  summaries, 
the  Manual  is  designed  to  serve  also  as  a  convenient  and  authoritative  reference  source. 
Those  already  familiar  with  a  topic  should  turn  immediately  to  the  summary  at  the  end  of 
the  appropriate  chapter.   They  can  determine  from  that  summary  what,  if  any,  of  the  body 
of  the  chapter  they  want  to  read  for  more  details. 

The  basic  approach  and  focal  point  of  the  treatment  in  this  Manual  is  the  measurement 
equation.   We  believe  that  every  measurement  problem  should  be  addressed  with  an  equation 
relating  the  quantity  desired  to  the  data  obtained  through  a  detailed  characterization  of 
the  instruments  used  and  the  radiation  field  observed,  in  terms  of  all  of  the  relevant 


parameters.   The  latter  always  include  the  radiation  parameters,  as  well  as  environmental 
and  instrumental  parameters,  as  previously  pointed  out.   The  objective  of  the  Manual  is  to 
develop  the  basic  concepts  and  characteristics  required  so  that  the  reader  will  be  able 
to  use  this  measurement-equation  approach.   It  is  our  belief  that  this  is  the  only  way  that 
uncertainties  in  the  measurement  of  optical  radiation  can  generally  be  limited  to  one,  or 
at  most  a  few,  per  cent. 

Currently,  the  Manual  deals  only  with  the  classical  radiometry  of  incoherent  radiation. 
The  basic  quantitative  relations  for  the  propagation  of  energy  by  coherent  radiation  (e.g., 
laser  beams)  are  just  being  worked  out  [19-22] . *   Without  that  basic  theory,  a  completely 
satisfactory  general  treatment  of  the  measurement  of  coherent  (including  partially  coherent) 
optical  radiation  is  not  possible.   Accordingly,  in  spite  of  the  urgent  need  for  improved 
measurements  of  laser  radiation,  we  won't  attempt  to  deal  with  it  now.   Possibly  this  sit- 
uation will  be  changed  before  the  current  effort  has  been  completed  and  a  supplement  on 
laser  measurements  can  be  added. 

As  stated  above,  we  first  hoped  to  prepare  this  Manual  on  a  more  elementary  level  but 
found  that  it  was  impossible  to  avoid  making  use  of  both  differential  and  integral  calculus 
of  more  than  one  variable.   However,  to  help  those  that  might  be  a  bit  "rusty"  with  such 
mathematics,  we  go  back  to  first  principles  each  time  a  mathematical  concept  or  procedure 
beyond  those  of  simple  algebra  or  trigonometry  is  introduced.   This  should  also  throw  addi- 
tional light  on  the  physical  and  geometrical  relationships  involved.   Where  it  seems  inap- 
propriate to  do  this  in  the  text,  we  cover  such  mathematical  considerations  in  appendices. 
It  is  also  assumed  that  the  reader  has  had  an  introductory  college  course  in  physics,  or 
the  equivalent. 

The  Manual  is  being  organized  into  three  Parts,  as  follows: 

Part  I.    Concepts 

Step  by  step  build  up  of  the  measurement  equation  in  terms  of  the  radiation  parame- 
ters, the  properties  and  characteristics  of  sources,  optical  paths,  and  receivers,  and  the 
environmental  and  instrumental  parameters.   Useful  quantities  are  defined  and  discussed 
and  their  relevance  to  various  applications  in  many  different  fields  (photometry,  heat- 
transfer  engineering,  astronomy,  photo-biology,  etc.)  is  indicated.   However,  discussions 
of  actual  devices  and  measurement  situations  in  this  Part  are  mainly  for  purposes  of 
illustrating  concepts  and  basic  principles. 

Part  II.    Instrumentation 

Descriptions,  properties,  and  other  pertinent  data  concerning  typical  instruments, 
devices,  and  components  involved  in  common  measurement  situations.   Included  is  material 
dealing  with  sources,  detectors,  filters,  atmospheric  paths,  choppers  (and  other  types  of 
optical  modulators),  prisms,  gratings,  polarizers,  radiometers,  photometers,  spectro- 
radiometers,  spectrophotometers,  etc. 

Part  III.    Applications 

Measurement  techniques  for  achieving  a  desired  level  of,  or  improving,  the  accuracy 
of  a  measurement.   Included  will  be  a  very  wide  variety  of  examples  of  environmental  and 
instrumental  parameters  with  discussion  of  their  effects  and  how  to  deal  with  them.   This 
is  where  we  deal  with  real  measurements  in  the  real  world.   The  examples  will  also  be 
drawn  from  the  widest  possible  variety  of  areas  of  application  in  illumination  engineering, 
radiative  heat  transfer,  military  infrared  devices,  remote  sensing,  meteorology,  astronomy, 
photo-chemistry  and  photo-biology,  etc. 


figures  in  brackets  indicate  literature  references  listed  at  the  end  of  this  Technical 
Note. 


Individual  chapter  headings  have  been  assigned  only  to  the  first  five  chapters: 

Chapter  1.    Introduction 

Chapter  2.    Distribution  of  Optical  Radiation  with  respect  to  Position  and 
Direction  —  Radiance 

Chapter  3.    Spectral  Distribution  of  Optical  Radiation 

Chapter  4.    Optical  Radiation  Measurements  —  a  Measurement  Equation 

Chapter  5.    More  on  the  Distribution  of  Optical  Radiation  with  respect  to 
Position  and  Direction 

Other  subjects  definitely  planned  for  Part  I  are  thermal  radiation,  photometry,  distribu- 
tion with  respect  to  time,  polarization,  diffraction,  and  detector  concepts.   It  is  not  our 
intention,  however,  to  try  to  complete  all  of  Part  I  before  going  on  to  Parts  II  and  III. 
In  fact,  because  we  realize  that  a  great  many  readers  are  probably  most  interested  in  the 
material  on  applications  to  appear  in  Part  III,  we  will  try  to  complete  and  publish  some 
chapters  in  Parts  II  and  III  just  as  soon  as  adequate  preparation  has  been  made  in  the 
earlier  chapters  of  Part  I.   However,  because  our  approach  to  radiometry  differs  so  much 
from  the  traditional  treatment,  we  feel  that  unnecessary  confusion  and  misunderstanding 
can  be  avoided  if  at  least  the  first  nine  chapters  of  Part  I  are  published  first  and  so 
are  available  to  readers  of  later  chapters. 

Finally,  we  invite  the  reader  to  submit  comments,  criticisms,  and  suggestions  for 
improving  future  chapters  in  this  Manual.   In  particular,  we  welcome  illustrative  examples 
and  problems  from  as  widely  different  areas  of  application  as  possible. 

As  previously  stated,  we  are  indebted  to  a  great  many  individuals  for  invaluable 
"feedback"  that  has  helped  us  to  put  this  text  together  more  effectively.   Notable  are  the 
inputs  and  encouragement  from  the  Council  on  Optical  Radiation  Measurements  (CORM) , 
especially  the  CORM  Coordinators,  Richard  J.  Becherer,  John  Eby,  Franc  Grum,  Alton  R. 
Karoli,  Edward  S.  Steeb,  and  Robert  B.  Watson,  and  the  Editor  of  Electro-Optical  Systems 
Design,   Robert  D.  Compton.   In  addition,  for  editorial  assistance,  we  are  grateful  to 
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SELF-STUDY  MANUAL  on  OPTICAL  RADIATION  MEASUREMENTS 
Part  I.    Concepts 

Chapter  6.    Distribution  of  Optical  Radiation  with  Respect  to  Polarization 

by  John  B.  Shumaker 

This  chapter  develops  and  illustrates  the  concepts  necessary 
to  include  polarization  rigorously  in  classical  radiometry.   The 
treatment  is  based  upon  the  Stokes  polarization  vector  of  spectral 
radiance  and  Mueller  transmittance  matrices.   The  Mueller  matrices 
of  many  common  polarizing  optical  components  are  discussed.   Con- 
siderable attention  is  paid  to  the  measurement  of  the  Stokes 
spectral-radiance-vector  components  and  to  the  measurement  of  the 
Mueller-transmittance-matrix  elements.   The  concepts  and  techniques 
are  illustrated  by  discussions  of  such  subjects  as  radiometer  cal- 
ibrations, three-polarizer  attenuators,  depolarizers,  and  the 
characterization  of  Polaroid-type  polarizers.   Also  included  are 
an  appendix  on  matrix  multiplication  and  an  appendix  showing  one 
way  by  which,  in  principle,  any  Mueller  transmittance  matrix  may 
be  measured. 

Key   Words:      Mueller  matrix;  polarization;  polarizer;  Stokes 
parameters. 

In  this  CHAPTER.    We  discuss  polarization,  another  of  the  radiation  parameters  introduced 
in  Chapter  1  [1].   The  other  radiation  parameters  —  position,  direction,  wavelength,  and 
time  —  are  continuous  variables  and  we  have  used  functional  notations  like  L(x,y,8 ,  <j>,  A) 
to  show  the  dependence  of  radiometric  quantities  upon  these  parameters.   Polarization,  how- 
ever, enters  somewhat  differently.   We  include  polarization  by  introducing  four  polariza- 
tion components  of  spectral  radiance,   L   n(x,y,8,<j>,  A)  ,   L   ,  (x,y,6,4>, A)  ,   L   „  (x,y,8  ,(j>,  A)  , 

A  ,  U  ^  ,  _1_  A  ,  Z 

and  L    (x,y,8,  <)>,  A) .   The  first  of  these,   L  n,      is  identical  with  the  spectral  radiance 
A  ,  j  A  ,  U 

we  have  been  using  in  earlier  chapters.   The  other  three,   L    ,   L    ,   and  L   Q,   pro- 

A,l     "  >  Z  "  ,  J 

vide  the  polarization  information.   The  expressions  that  result  from  this  four-fold 
extension  of  the  quantities  appearing  in  the  measurement  equation  are  adequate  to  handle 
any  kind  of  polarizing  optical  element  and  any  polarization  state  within  the  approximations 
of  classical  (geometrical  optics)  radiometry  of  incoherent  radiation. 

We  will  not  in  this  chapter  discuss  in  detail  specific  polarizing  or  depolarizing 
instruments.   These  will  be  described  in  Part  II  and  are  also  covered  in  several  of  the 
references  [2,3,4,5].   We  concentrate  here  on  developing  and  illustrating  the  use  of  the 
concepts  and  equations  necessary  to  include  polarization  rigorously  in  classical  radiom- 
etry.  We  begin  the  chapter  with  a  review  of  the  relationship  between  the  oscillating 
electric  field  and  the  polarization  state  of  radiation.   We  then  introduce  the  four  com- 
ponents of  spectral  radiance  and  study  how  they  are  transformed  by  ideal  polarizers  and 
retarders.   Next  we  discuss  the  response  of  a  radiometer  to  polarized  light  and  how  it  may 
be  calibrated.   After  this  we  examine  real,  non-ideal,  optical  components  and  finally  con- 
clude the  chapter  with  a  brief  look  at  some  experimental  precautions  in  the  use  of  polar- 
izing optics. 


The  object  of  this  chapter  is  to  extend  the  equations  of  radiometry  which  have  been 
developed  in  preceding  chapters  so  that  they  can  include  polarization  rigorously.   We  make 
little  attempt  to  provide  an  intuitive  understanding  of  the  polarization  phenomenon.   Our 
approach  is  directed  toward  the  radiometrist  who  may  regard  polarization  as  an  unfortunate 
complication  which,  like  the  common  cold,  has  been  inflicted  upon  him  by  the  whimsical 
gods.   We  therefore  have  tried  to  present  a  systematic  measurement  approach  to  the  polari- 
zation problem  which  does  not  require  the  reader  to  become  an  expert  on  polarization  optics. 
As  a  consequence,  we  illustrate  our  equations  and  the  kinds  of  measurement  information 
involved  with  conceptual  measurements  which  employ  simple,  easily  available  instruments. 
Although  some  very  precise  relative  polarization  measurements  can  be  made  using  high 
quality  optical  components  and  specialized  instruments  such  as  ellipsometers,  such  preci- 
sion is  usually  not  required  in  radiometry.   These  special  devices  and  techniques  are 
described  in  some  of  the  references  [2,5,6]  and  will  also  be  considered  in  Parts  II  and  III. 
This  chapter  also  does  not  discuss  the  many  remarkable  uses  which  man  and  nature  have  made 
of  polarized  light.   These  are  briefly  discussed  in  reference  [2]  which  includes  extensive 
further  literature  references. 

THE  STOKES  COMPONENTS  of  POLARIZATION 

Polarization  can  be  thought  of  in  terms  of  the  orientation  of  an  electromagnetic  field 
vector  oscillating  about  the  line  of  sight  between  an  observer  and  a  light  source.   All  of 
the  consequences  of  polarization  can  be  derived  from  this  starting  point  and  the  names 
"linear"  polarization  and  "circular"  polarization  originate  in  this  description.   However, 
within  the  approximations  of  geometrical  optics  a  more  direct  phenomenological  approach  is 
possible.   This  simply  requires  the  reinterpretation  of  spectral  radiance  as  being  fully 
characterized  by  four  independent  polarization  components,  called  Stokes  components.   All 
four  components  are  measurable,  as  we  will  see,  and  form  a  phenomenological  basis  for  the 
quantitative  inclusion  of  polarization  in  radiometry  and  photometry.   There  are  other  ways 
of  describing  polarization  [2,4,5]  but  for  our  purposes  the  use  of  Stokes  components  is  the 
simplest  because  for  light  with  a  negligible  degree  of  polarization  —  frequently  the  most 
important  situation  in  radiometry  —  the  treatment  reduces  directly  to  the  equations  we 
have  been  considering  in  the  preceding  chapters.   The  Stokes  components  provide  no  particu- 
lar insight  into  why  or  how  polarizing  optical  components  work  but  we  will  not  be  inter- 
ested in  that  aspect  of  polarization  in  this  chapter.   We  are  interested  in  how  the  polar- 
ization state  of  a  beam  affects  the  response  of  a  radiometer  and  in  how  to  describe  the 
flux  of  a  polarized  beam  so  that  the  beam  is  uniquely  defined  from  a  radiometric  point  of 
view  —  so  that  it  could  be  used  to  calibrate  a  radiometer,  for  example.   The  four  Stokes 
components  of  spectral  radiance  at  every  point  and  direction  in  a  radiation  field  provide 
that  kind  of  description.   Since  most  people  prefer  to  think  in  the  familiar  terms  of 
electric  field  vectors  and  linear  and  circular  polarization,  we  will  indicate  how  the 
Stokes  components  are  related  to  these  common  concepts,  but  we  emphasize  that  this  is  not 
essential  to  radiometry. 


We  begin  with  the  familiar  microscopic  description  of  polarization.   We  consider  a 
single  steadily  radiating  oscillator  at  a  distant  point  on  the  negative  Z-axis  of  a  rectan- 
gular coordinate  system  and  imagine  that  we  can  continuously  record  the  instantaneous 
electric  field1  present  at  the  origin  due  to  this  oscillator.   Except  at  distances  from  the 
light  source  of  a  few  wavelengths  or  less,  the  Z  component  of  the  electric  field  —  the 
component  along  the  direction  of  propagation  —  will  be  completely  negligible  and  the  field 
will  lie  in  the  X-Y  plane.   The  X  and  Y  field  components  then  will  be  of  the  form 

E     =  V    'cosUir'vt  +   6    ) 

X  X  X  [V-nT1]  (6.1) 

E     =  V    -cos(27T-vt  +  6    ) 

y        y  y 

where  V   and  V   are  the  amplitudes   [V'm-1],   v   is  the  frequency   [Hz],   6   and   6 
x        y  x        y 

are  the  phases  trad]  of  the  electromagnetic  wave  train,  and   t   is  the  time   [s].   In 

figure  6.1  are  shown  some  examples  to  illustrate  how  the  electric-field  vector  changes  with 

time  for  some  special  combinations  of  amplitudes  and  phases.   In  each  case  the  sketch  on 

the  left  shows  the  pattern  traced  out  in  the  X-Y  plane  by  the  tip  of  the  electric-field 

vector  and  the  graphs  on  the  right  show  the  time  behavior  of  the  X  and  Y  components.   In 

diagrams  (a)  we  see  light  which  is  described  as  linearly  polarized2  in  a  vertical  direction; 

the  horizontal  component,   E  ,   of  the  electric  field  is  always  zero.    Diagrams  (b)  show  a 

more  general  example  of  linear  polarization.   E   and  E   are  exactly  in  phase  with  one 

another,   6   =  6  ,   but  their  amplitudes  V   and  V   differ.   Diagrams  (c)  show  right 
y    x  y        x  — ° — 

circularly  polarized  light;  the  electric  vector  appears  to  trace  a  circle  in  the  clockwise 


The  electric  and  magnetic  field  vectors  are  perpendicular  to  each  other  and  to  the  direc- 
tion of  propagation  of  the  radiation.   Our  development  could  equally  well  have  been  based 
upon  the  magnetic  field,  as  it  is  in  much  of  the  older  literature. 

Linear"  and  "plane"  polarization  are  synonymous  terms.   "Linear"  seems  more  consistent 
with  the  terms  "circular"  and  "elliptical"  and  so  is  used  exclusively  here. 

3The  terms  "horizontal"  and  "vertical"  linear  polarization  and  "right"  and  "left"  circular 
polarization  are  somewhat  arbitrary  and  are  used  differently  by  different  authors.   Light 
which  contains  no  vertical  component  of  the  electric  field  may  be  described  as  horizontally 
polarized;  however,  if  the  magnetic  field  vector  is  emphasized  instead  of  the  electric 
vector  then  such  light  would  be  called  vertically  polarized.   Right  and  left  handedness  are 
likewise  ambiguous  descriptions  of  circular  polarization.   Fortunately,  as  long  as  one  is 
consistent  within  his  own  calculations  it  rarely  matters  what  conventions  one  adopts  for 
these  definitions.   To  insure  consistency  one  must  usually  take  all  his  formulas  and  equa- 
tions from  the  same  source.   Our  conventions  used  in  this  chapter  agree  with  those  of 
Shurcliff  [2]  where,  for  instance,  light  reflected  from  horizontal  surfaces  such  as  from  a 
desk  top  or  the  ocean  is  defined  as  possessing  enhanced  horizontal  polarization.   Polarizing 
sunglasses,  consequently,  are  vertical  polarizers  as  customarily  worn. 


Figure  6.1 


Electric-field-vector  behavior  for 
monochromatic  polarized  radiation 


(a)   Linearly  polarized  in  the  Y-Z  plane 


(b)   Linearly  polarized  at  +30' 


-^v 


(c)   Right  circularly  polarized 


(d)   Right  elliptically  polarized 


direction  (as  we  look  back  toward  the  source) .   E   and   E   have  equal  amplitudes  but  are 

y        x 

exactly  a  quarter  cycle  out  of  phase  with  one  another:   6=6+  it/2.   If   6   =  6   -  tt/2 

y    x  y    x 

the  electric  vector  will  rotate  in  a  counterclockwise  direction  and  the  polarization  is 
called  left  circular.   The  final  example  shows  elliptically  polarized  light  which  is  the 
general  case  of  eq.  (6.1).   In  this  case  no  special  relationships  exist  between  V   and 
V   and  between   6    and   6  .   E   and   E   have  unequal  amplitudes  and  are  out  of  phase 
by  an  angle  y  =  6   -  6  .    If  y      is  positive  the  polarization  is  called  right  elliptical 

y   x 

and  if  y      is  negative  it  is  called  left  elliptical.   Elliptical  polarization  represents 
the  most  general  expression  of  a  perfectly  monochromatic  ray  which  is  possible.   Linear  and 
circular  polarizations  are  special  cases  of  elliptical  polarization. 

As  an  aid  to  visualizing  polarization  we  might  think  of  a  classical  source  mechanism 
in  terms  of  the  oscillations  of  charged  particles.   At  radio  wavelengths  a  vertically 
polarized  field  can  be  produced  by  a  vertical  transmitting  antenna  with  a  cluster  of  elec- 
trons sloshing  up  and  down  at  the  radio  frequency.   At  optical  wavelengths  we  can  similarly 
imagine  a  vertical  vibration  of  electrons  at  the  source  as  one  way  of  generating  vertically 
polarized  light.   Circular  polarization  then  can  be  thought  of  as  being  generated  by 
charges  moving  in  circles  at  the  source.   If  the  charges  are  moving  in  circular  orbits  as 
seen  by  an  observer  in  one  direction  then  observers  in  other  directions  will  see  the 
charges  moving  in  ellipses  or  even  vibrating  in  straight  lines  thus  observing  elliptical  or 
linear  polarization.   Such  directional  dependence  of  the  polarization  form  is  actually 
observed  in  the  Zeeman  effect,  for  example,  where  some  spectral  lines  emitted  by  a  gas  in  a 
magnetic  field  show  circular  polarization  if  viewed  along  the  direction  of  the  magnetic 
field  and  change  continuously  through  elliptical  to  linear  polarization  as  the  viewing 
direction  changes  toward  the  perpendicular  to  the  magnetic  field. 

Three  parameters  are  required  to  describe  the  general  polarization  ellipse.   They  can 
be  chosen  in  several  ways.   For  example  the  magnitude  of  the  parameter  y  =   6   -  6 

y   x 

specifies  the  ellipticity  of  the  polarization  and  its  sign  specifies  its  handedness  —  left 

or  right.   The  parameter  V /V   specifies  the  orientation  of  the  ellipse  in  the  X-Y  plane. 

And  the  parameter  V  2  +  V  2   specifies  the  size  of  the  ellipse.   A  set  of  parameters  which 
x     y 

has  proven  to  be  somewhat  more  useful,  however,  is  the  set  of  Stokes  parameters2: 


Y   is  further  required  to  be  adjusted  by  adding  or  subtracting  integral  multiples  of   2tt 
if  necessary  to  give   |y|  <  ^ • 

In  eq.  (6.2)  that  follows,  the  proportionality  constant   k   is  given  by  k  =   n/(2yc) 
where  n   is  the  refractive  index,   u   the  magnetic  permeability  and   c   is  the  vacuum 
velocity  of  light.   If  V   and  V   are  expressed  in  volts/meter  and   E   *   in  watts/m2 
then  </n  =  1.327  x  10~3  [W*V-2]  for  non-magnetic  materials.   Equations  (6.2)  contain  less 
information  than  eqs.  (6.1).   The  absolute  phase  information  has  been  lost.   Ultimately 
this  will  result  in  the  exclusion  of  interference  effects  from  our  treatment  and  in  the 
additivity  of  radiometric  quantities  such  as  spectral  radiance.   These,  of  course,  are 
just  the  limitations  (and  advantages)  of  geometrical  optics. 

5 


E  n   =  k(V  2  +  V  2)    [W-m"2] 
n,0      x     y 

E   .  =  k(V  2  -  V  2) 
n,l      x     y 

E   .  =  k(2V  -V  -cosy) 
n,2       x  y 

E   _  =  k(2V  -V  -siny). 
n,3       x  y 


(6.2) 


There  are  four  of  these  parameters,  the  first  of  which  is  redundant  in  the  present  context1 

but,  as  we  shall  see  later,  forms  the  basis  for  including  unpolarized  light.   The  second 

Stokes  parameter  E  -.      is  a  measure  of  the  excess  of  horizontal  linear  polarization  over 
n,l 

vertical  linear  polarization.   It  is  called  the  horizontal  preference  parameter.   The  third 

parameter   E   .   is  a  measure  of  the  excess  of  linear  polarization  at  45°  (from  the  hori- 
n,  i. 

zontal)  to  that  at  135°.   It  is  called  the  +45°  preference  parameter.   The  last  Stokes 
parameter  is  a  measure  of  the  excess  of  right-circular  component  to  left-circular  component. 
It  is  called  the  right-circular  preference  parameter.   These  parameters  all  have  the  unit- 
dimensions  of  irradiance  and  constitute  the  four  Stokes  polarization  components  of  normal 
irradiance  from  a  single  steady  oscillator. 

A  real,  non-monochromatic,  non-point  radiation  source  will  also  generate  a  varying 
electric  field  vector  at  any  point.   But  the  motion  of  this  electric  field  vector  will 
generally  not  be  easy  to  describe.   For  one  thing  the  motion  need  not  be  confined  to  a 
plane.   For  another,  it  may  not  be  periodic.   In  the  limit  as  the  solid  angle  and  wavelength 
interval  approach  zero  the  motion  of  the  electric  field  vector  will  become  planar  and 
elliptical  and  so  the  Stokes  parameters  are  appropriate  for  describing  the  polarization 
state  of  the  radiometric  quantity  spectral  radiance.   Accordingly,  the  average  values  of  the 
Stokes  spectral  radiance  components  that  a  real  radiometer  would  measure  are  given  by 
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(6.3) 


A,  3   AA-At 


AAJ 


At 


2V  -V  -siny 
x  y    ' 

Aoo-AA 


dx-dy-dt 


where  At   is  the  time  required  for  a  measurement,   AA  is  the  area  of  the  detector  surface 
perpendicular  to  the  ray  direction,   Aw   is  the  solid  angle  from  which  rays  can  illuminate 


because     E2=E2+E2+E2 
n,0  n,l  n,2  n,3 


the  detector  and   AA   is  the  wavelength  band  reaching  the  detector.   V  ,   V   and   y   are 

x    y 

the  field  amplitudes  and  phase  difference  at  the  detector  surface. 

In  the  mathematical  limit  as   AA,   At,   Aw  and   AA   approach  zero 

.       LA2,0  "   LA2,1  +  LA2,2  +  h2  ,3-  (6'4) 

This  is  the  condition  for  complete  polarization.   With  sufficient  radiometer  resolution  in 
these  four  parameters,   A,   t,   oj  and   A,   all  measurement  results  would  show  complete 
polarization.   However,  as  was  discussed  in  Chapter  2,  [1]  radiometric  measurements  cannot 
be  made  if  any  of  these  resolution  intervals  goes  to  zero  because  zero  power  would  then  be 
received.    The  practical  limits  on  the  smallness  of   AA,   At,   Aid   and   AA   result  in 
values  of  Stokes  spectral-radiance  components  which  obey 

LA2,0  *  LA2,1  +  LA2,2  +  LA2,3' 
Thus,  we  speak  of  a  polarized  component  of  spectral  radiance 


and  an  unpolarized  component 


LA,P   =    VLA2,1  +  L^2  +  LA2,3  (6'5) 


L  =   L  -   L. 

A  ,u  A  ,0  A,p 


and   define   the  degree   of   polarization  as 

P   "   LA,P/LA,0-  (6-6) 

Although  at  first  glance  it  would  seem  that  a  definition  of  unpolarized  light  which 
obviously  depends  upon  instrument  resolution  would  be  highly  unsatisfactory,  it  turns  out 
that  this  is  rarely  a  problem.   Usually  unpolarized  light  arises  because  of  inadequate 


By  considering  the  coherence  properties  of  radiation  one  can  deduce  that  the  state  of 
polarization  in  the  detector  plane  remains  essentially  constant  over  areas  of  the  order  of 
6A  'v  A2*  (6oj)-1  .   Similarly,  the  state  of  polarization  is  essentially  constant  in  time  over 
a  time  interval  of  the  order  of   6t  ^  (6v)_1  =  A2*(c*6A)_1.   Therefore,  the  resolution 
intervals  need  be  no  smaller  than  these  values.   A  perfect  radiometer  with  these  resolution 
intervals  would  require  a  source  with  spectral  radiance 


h*c    / , 9    A 


L,    n  *   104   h-v/(6u)-6A-6A-5t)    =   lO4  •—.—  •(  A2-—)        <\.   2-1019         [W-m-3 -sr"1 ] 
A  ,1)  Aye/ 

to  achieve  a  1%  measurement  accuracy  at  A  =  500  [nm] .  This  is  a  factor  of  106  times  the 
solar  spectral  radiance.  For  any  less  intense  source  the  radiometer  resolution  intervals 
must  be  larger  if  useful  accuracy  is  desired. 


resolution  in  the  time  dimension.   For  typical  spectrometers  and  for  steady  incandescent 

sources  (so-called  natural  light),  measurement  times  less  than  10-11  seconds  would  be 

required  to  observe  complete  polarization.   That  is,   V  ,   V  ,   and   y  will  not  change 

x    y 

appreciably  over  periods  of  10-11  seconds  and  non-vanishing  values  of  L ,  ,,   L   „   or 

A ,  1     A  ,  z 

L,  „   can  be  expected.   If   At   is  very  much  longer,  however,  the  fluctuations  in  V  ,   V 
A, J  x    y 

and   y  may  cause  cancellations  in  computing  the  mean  values  of  these  three  Stokes  compo- 
nents.  L,  „,   being  the  sum  of  positive  quantities   V2  and  V2,   will  not  be  reduced  by 
A,0  x        y 

averaging.   Thus,   L     will  appear  to  include  a  contribution  which  has  disappeared  from 
A ,  U 

the  L    ,   L   „   and  L     components  and  is  consequently  called  unpolarized.   Since  the 

A,i     A,z  A, j 

time  resolution  required  to  resolve  the  unpolarized  component  is  so  far  beyond  the  capa- 
bilities of  customary  radiometry  there  is,  for  all  practical  purposes,  no  ambiguity  in  the 
definition  of  unpolarized  light.   Occasionally  a  (usually  periodic)  rapid  variation  in 
polarization  over   t   or  over  one  of  the  other  parameters  A,   oj  or   A   is  deliberately 
introduced  to  simulate  the  "natural"  unpolarized  light  we  have  just  described.   These 
variations  are  usually  sufficiently  coarse  that  they  could  conceivably  be  resolved  by  some 
radiometers.   For  this  reason  such  light  is  frequently  called  depolarized  instead  of 
unpolarized. 

The  first  Stokes  component,   L  n,   is  the  spectral  radiance  that  we  have  been  con- 

A ,  u 

sidering  in  earlier  chapters.   It  is  a  measure  of  the  total  radiant  flux  associated  with 

a  ray  and  is  therefore  usually  the  quantity  of  greatest  interest.   It  is  the  value  of 

spectral  radiance  which  a  polarization-indifferent  radiometer  would  measure.   L  n   can 

A  ,  U 

never  be  negative.   The  other  three  Stokes  components   L ,  , ,   L   „   and  L     may  be 

A, I     A  ,  z         A ,  j 

either  negative  or  positive.   If  the  radiation  contains  more  vertical  linear  polarization 

than  horizontal  linear  polarization  then  the  horizontal  preference  component  L     will  be 

A, I 

negative.   If  the  radiation  is  completely  vertically  polarized  then  L   ..  =  -L    .   If  the 

A  ,  I      A  ,  u 

radiation  contains  equal  horizontally  and  vertically  linearly  polarized  contributions  as  in 

unpolarized  light  or  in  45°  polarized  light  then  L    =  0.   If  the  horizontal  content  is 

A,l 

in  excess  then  L     is  positive;  its  maximum  value  is   L     when  the  radiation  is  com- 
A  ,  1  A  ,  U 

pletely  horizontally  linearly  polarized.   The  +45°-preference  Stokes  component  of  spectral 

radiance,   L   „,   describes  the  +45°   and  +135°   linearly  polarized  content  of  the  radia- 

A  ,z 

tion  in  an  exactly  analogous  manner.   And  finally,  the  right-circular  preference  component, 

L    ,   similarly  takes  on  values  from  -L   _   for  completely  lef t-circularly  polarized 
A  ,  3  A  ,  U 

light  through  zero  for  light  with  no  net  circular  polarization  to  +L   _   for  completely 

A  ,  U 

right-circularly  polarized  light.   Radiation  linearly  polarized  at  other  angles  (other  than 
0°   and   45°   and  the  complementary  angles   90°   and  135°)   will  have  non-zero  values  for 
both  L,  .   and   L.  „.   Likewise  if  either  or  both  of   L,  ..   and  L   0   are  non-zero  and  in 

A, I  A,Z  A ,  J.  A,Z 

addition  L     is  non-zero  then  the  radiation  contains  attributes  of  both  linear  and  cir- 

A,  3 

cular  polarization  and  is,  of  course,  elliptically  polarized. 

Although  for  spectral  radiance  the  four  Stokes  components  provide  a  complete  and,  for 
all  practical  purposes,  unique  description  of  the  polarization  state  of  radiation,  the  same 
cannot  be  said  for  the  integrals  of  spectral  radiance  such  as  radiance,  spectral  intensity, 
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spectral  irradiance,  etc.   For  such  cases,  as  we  have  indicated,  the  trajectory  of  the  tip 
of  the  electric  vector  can  no  longer  be  described,  in  general,  by  three  parameters.   The 
only  Stokes  component  of  these  integrated  radiometric  quantities  which  is  always  unambigu- 
ous is  the  first  one  —  the  integrals  of   L   „.   As  in  earlier  chapters  where  polarization 

A  ,  U 

was  neglected,  an  integral  of  spectral  radiance  does  not  uniquely  characterize  a  complete 

radiation  field;  one  must,  in  general,  trace  the  spectral  radiance  by  determining  the 

propagance,  ray  by  ray,  from  source  to  receiver  and  then,  at  the  receiver,  perform  the 

desired  integration  of  spectral  radiance.   Fortunately,  at  the  ultimate  detector  the  usual 

quantity  of  interest  is  the  total  power,  and  to  calculate  this  requires  only  an  appropriate 

integral  of   L    .   Therefore,  we  shall  discuss  polarization  only  in  terms  of  rays  and 
A  ,  U 

spectral  radiance. 

STOKES  VECTORS  and  MUELLER  MATRICES 

The  complete  description  of  polarization  in  classical  radiometry  requires  the  specifi- 
cation of  four  parameters  such  as  the  Stokes  components  of  spectral  radiance  associated  with 
each  ray.   These  Stokes  components  are  measurable  with,  for  example,  a  linear  polarizer  and 
a  quarter -wave  plate.   We  will  return  to  this  later  after  we  have  generalized  the  concept 
of  propagance  or  transmittance  to  include  polarization. 

In  Chapter  3,  eq.  (3.14)  [1]  we  related2  the  spectral  radiance  L'   of  a  ray  emerging 
from  an  optical  medium  to  the  incident  spectral  radiance  L  by  means  of  a  transmittance  or 
propagance  T: 


^For  simplicity  we'll  usually  speak  only  of  transmittance  and  won't  distinguish  between 
propagance  and  transmittance.   However,  our  results  apply  equally  well  when  reflectance  is 
involved  provided  some  care  is  exercised  in  defining  the  coordinate  system  used  for  the 
incident  and  reflected  rays.   The  simplest  coordinate  system  to  use  is  one  with  the  origin 
at  the  point  of  reflection  of  the  ray.   The  incident  ray  is  referred  to  such  a  coordinate 
system  with  its  Z-axis  along  the  ray  in  the  direction  of  propagation,  its  Y-axis  in  the 
plane  of  incidence  —  that  is,  the  plane  defined  by  the  incident  and  reflected  rays  —  and 
its  X-axis  perpendicular  to  the  plane  of  incidence  and,  therefore,  lying  in  the  reflecting 
surface  or,  more  generally,  tangent  to  it.   The  reflected  ray  is  referred  to  a  coordinate 
system  with  its  Z-axis  in  the  direction  of  propagation  of  the  reflected  ray,  with* its  X-axis 
coincident  with  the  X-axis  of  the  incident  ray  and  with  its  Y-axis  in  the  plane  of  incidence 
but  with  its  sense  chosen  to  preserve  the  handedness  of  the  coordinate  system.   If  the 
incident-ray  coordinate  system  is  right-handed  then  that  of  the  reflected  ray  is  also. 

2In  the  interests  of  typographical  simplicity  several  of  the  symbols  used  in  earlier  chap- 
ters will  be  modified  in  this  chapter.   In  particular,  since  we  will  be  dealing  only  with 
spectral  radiance  we  will  have  no  need  to  distinguish  between  L   and  L  ,   so  we  will 

A 

usually  omit  the  subscript  A  on  this  quantity.  The  spectral-radiance  input  to  a  device, 
formerly  denoted  L  . ,  will  usually  be  written  L  and  the  output,  L  ,  will  usually  be 
written  L' .   The  transmittance,   t,   or  propagance,   t* ,   will  usually  be  written  T. 
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L'(x,y,e,<M)  =  T(x,y,0,(f>,A)-L(x,y,e,(t>,A). 

L  and  L'   we  now  see  consist  of  four  components  each  so  we  must  expect  four  equations  of 
this  type  —  one  for  each  component.   Moreover,  since  the  propagance  may  differ  for  differ- 
ent forms  of  polarization  and  since  every  component  of  L  may  in  general  contribute  to 
each  component  of   L'   the  propagance   T  must  now  contain  sixteen  elements: 

L0  =  T00'L0  +  T0rLl  +  T02*L2  +  T03'L3 


Ll  =  T10*L0  +  Tll'Ll  +  T12"L2  +  T13'L3 


L2  "  T20-L0  +  T21'L1  +  T22*L2  +  T23'L3 


L3     T30*L0  +  T31*L1  +  T32*L2  +  T33*L3" 


(6.7) 


This  is  the  generalization  of  eq.  (3.14)  which  is  required  to  describe  the  changes  in  the 
Stokes  spectral-radiance  components  of  a  ray  upon  passing  through  an  optical  medium.   Just 
as  the  spectral  radiance  in  a  beam  may  vary  from  ray  to  ray  and  from  wavelength  to  wave- 
length, so  may  the  propagance  elements   Tnn,   Tm ,  .  .  .  vary  from  ray  to  ray,  etc.   Thus 
the  values  of  these  elements  depend  upon  the  ray  we  are  examining  and  we  should  write  them 
as   Tnn(x,y,  6,  <j>,  A) ,   etc.,  where  the  x,y,6,   and   <j>  are  the  same  reference-surface 

quantities  which  identify  the  ray  of  spectral  radiance  L   n(x,y,6  ,cj>,  A)  .   However,  through- 

A  ,  U 

out  most  of  this  chapter  we  deal  with  the  behavior  of  a  single  ray  so  we  will  usually  not 
explicitly  display  this  dependence. 


Because  the  four  Stokes  polarization  components  are  independent  [within  the  restric- 
LQ>  (L^  +L^  +L^] 
Stokes  spectral-radiance  vector 


tion  L   >  (L2  +  L  2  +  L  2)  2  ]   it  is  common  to  describe  the  set  of  four  as  a  vector,  the 


■L   = 


~Lo 

Ll 

L2 

_L3 

(6.8) 


By  expressing  spectral  radiance  this  way,  as  a  column  vector,  eqs.  (6.7)  can  be  compactly 
written  as 


L'    -  T  •  L  . 


(6.9) 


where  is   the      4x4  matrix 


10 


T    = 


T   T    T    T 
00   01   02   03 

T   T    T    T 
10   11   12   13 

T    T    T    T 
20  21  22  23 

T   T   T    T 
L  30  31  32   33 


(6.10) 


and  eq.  (6.9)  means  exactly  the  same  thing  as  eqs.  (6.7)  —  no  more  and  no  less.   We  shall 
generally  use  this  vector  and  matrix  notation1  throughout  the  rest  of  this  chapter,  not 
because  of  any  useful  four-dimensional  mental  image  that  might  be  envisioned,  but  simply  for 
economy  of  notation. 

There  is  one  precaution  about  matrix  equations  which  must  be  observed.   That  is  that 
the  product   f  •  |_  is  not  always  the  same  as    •  ]   and  the  product  of  two  matrices 
Ti  "To   i-s  not  generally  the  same  as  To*T -i  •   Tne  rule  for  writing  the  equations  correctly 
in  accordance  with  the  conventional  matrix  multiplication  formulas  is  to  write  the  incident 
radiation  vector  and  the  transmittance  matrices  in  the  order,  from  right  to  left,  in  which 
the  ray  encounters  the  optical  elements  they  represent.   So,  if  a  ray  with  spectral- 
radiance  vector  |_  passes  through  one  optical  device  with  transmittance  matrix  J.,   and 
then  through  a  second  with  transmittance  matrix  J„ ,   the  emerging  spectral  radiance  will 
be 


L'(x,y,9,(j),A)  =  T2(x,y,6,<»,A)  •  ^(x.y.e.^.A)  •  |_(x,y,e,<f>,A) 


(6.11) 


The  effect  of  a  succession  of  optical  devices  can  be  combined  into  a  single  effective  trans- 
mittance matrix  by  simply  multiplying  the  individual  matrices  together  in  the  proper  order. 
Thus,  eq.  (6.11)  can  be  written 


L'    =   T-L  . 


(6.12) 


where  the  product  matrix   f  =  To'Tt   can  t>e  obtained  by  following  the  rules  for  matrix 
multiplication.   Of  course  the  same  result  would  be  obtained  by  writing  out  eqs.  (6.7)  for 
the  radiance  components  passing  device  1  and  then  again  using  eqs.  (6.7)  inserting  these 
expressions  for  the  radiance  input  to  device  2  to  obtain  expressions  for  the  final  emerging 
radiance  components. 

The  matrices  which  describe  the  transmittance  or  propagance  through  optical  elements 
which  influence  polarization  are  known  in  the  literature  as  Mueller  matrices.   In  very 
simple  cases  the  proper  form  of  the  Mueller  matrix  is  obvious  by  inspection.   For  optical 
elements  for  which  equations  can  be  written  relating  the  amplitudes  and  phases  of  the 


headers  who  are  not  familiar  with  matrix  manipulations,  or  who  tend  to  forget  which  sub- 
script means  what,  will  find  everything  needed  to  convert  these  matrix  equations  to  workable 
formulas  in  Appendix  6. 
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emerging  radiation  field  to  the  amplitudes  and  phases  of  the  entering  field,  the  matrix  can 
be  derived  by  a  straightforward  procedure  [7,8].   This  includes  all  perfect  polarizing 
optical  elements  such  as  ideal  polarizers,  quarter-wave  plates,  mirrors,  etc.   For  a  real, 
non-ideal  device,  however,  the  Mueller  matrix  can  only  be  obtained  ultimately  by  experiment. 
Mueller  matrices  for  ideal  optical  components  have  been  tabulated  by  many  authors  [2,7]. 

THE  SIMPLE  SPECTRAL  FILTER 

Let  us  now  study  some  examples  of  the  use  of  Mueller  matrices  and  Stokes  components. 
Perhaps  the  simplest  useful  optical  element  is  just  a  glass  absorbing  spectral  filter  with 
transmittance   t(A).   The  Mueller  transmittance  matrix  for  this  should  have  the  effect  of 
simply  multiplying  each  component  of   |_  by   t   with  no  mixing  by  or  interaction  with  other 
components.   That  is,  we  want 


Lo' 

=   x(A)-L 

Ll 

-   t(A)-L 

L2 

=   x(A)-L 

L3 

=   t(A)-L 

(6.13) 


Comparing  these  equations  with  eq.  (6.7)  we  see  that   T 


00 


x(A),   Tn  =  t(A), 


T„_  =  t(A)   and  all  other   T   =  0.   Thus,  for  a  simple  filter  we  have 


T22=  x(A), 


T   = 


t(A)  0  0  0 

0  t(A)  0  0 

0  0  t(A)  0 

0  0  0  t(A) 


=  t(A) 


1 

0 

0 

0 

0 

1 

0 

0 

0 

0 

1 

0 

0 

0 

0 

1 

(6.14) 


THE  IDEAL  POLARIZER 

A  more  interesting  optical  element  is  the  ideal  linear  polarizer.   An  ideal  polarizer 
possesses  a  polarization  axis  (perpendicular  to  the  ray  direction)  with  the  property  that 
the  transmittance  to  light  linearly  polarized  in  this  direction  is  1  and  the  transmittance 
to  light  linearly  polarized  at  right  angles  to  this  direction  is  0.   Aside  from  surface- 
reflection  effects,  crystal-prism  polarizers  such  as  the  Nicol-prism  polarizer  come  close 
to  this  ideal.   If  the  angle  between  the  polarization  axis  of  the  polarizer  and  the  hori- 
zontal reference  axis  or  x-axis  of  the  optical  system  is   <j>   the  Mueller  transmittance 
matrix  of  an  ideal  linear  polarizer  is 


POD)     = 


1  cos2tf>  sin2<f>          0 

cos2tf)  cos22<j)  sin2<f>*cos2<j)   0 

sin2<|>  sin2(f>' cos 2(f)  sin22cj>          0 

0  0  0               0 


(6.15) 
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Let  us  use  eq.  (6.9)  to  see  the  effect  of  such  a  polarizer  on  unpolarized  light.   Since 
unpolarized  light  is  characterized  by  L^  =  L2  =  L3  =  0,   we  have 


L'  =  P(*)-L  =  h 


1 

cos2<J> 

sin2<J> 

0 

L 

cos2$ 

cos   2<j> 

sin2<f>'cos2tj> 

0 

0 

sin2<fi 

sin2(j)*cos2(}) 

sin22<f> 

0 

0 

0 

0 

0 

0 

0 

0 


=  h  l. 


1 

cos2c| 

sin2<| 

0 


(6.16) 


The  vector   equation    (6.16)    is   simply  a   shorthand   for   the   set   of   four   equations: 


Lx'     =  JiL0-cos2<j) 
L2'     =  %L0-sin2<j> 

L3     "  °- 

By  using  the  trigonometric  identity   cos22<j)  +  sin22cf)  =  1  we  see  that 

L'  =  (l'2  +  L'2  +  L'2)1,   This  shows  [eq.  (6.4)]  that  the  output  radiance  is  completely 


polarized.   Since  L  '  =  0   there  is  no  circular  component  so  the  polarization  is  strictly 
linear.   If  we  turn  the  polarizer  so  that  its  axis  is  parallel  to  the  x  axis  we  would 
expect  horizontally  polarized  light.   Putting   <j>  =  0°   in  eq.  (6.16)  we  obtain 


L'  =  P(o°)-L  =  k\ 


Thus,  light  with  Stokes  components  like 


is  pure  horizontally  polarized  light.   If  the  polarizer  is  turned  to 
light  polarized  in  that  orientation.   From  eq.  (6.16)  we  get 


=  45°  we  expect 


L'  =  hLQ 


Thus,  light  whose  Stokes  components  are  a  multiple  of 
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is  linearly  polarized  at  45°.   Similarly  inserting   if)  =  90°   and   135°   in  eq.  (6.16) 
shows  that  vertically  polarized  light  is  a  multiple  of 


1 

-1 

0 

0 


and  light  polarized  at  135°  is  characterized  by  Stokes  components  proportional  to 


1 

0 

-1 

0 


Notice  that  the  spectral  radiance  L'    is  always  positive,  as  its  name  suggests  it  must  be, 
while  the  other  components  may  have  either  sign. 

Suppose  instead  of  starting  with  unpolarized  light  we  had  put  arbitrarily  polarized 
light  through  the  ideal  linear  polarizer.   We  have 


L'  =  P(*)-L  =  h 


1 

c 

s 

0 

Lo 

c 

c2 

s-c 

0 

Ll 

s 

s-c 

s2 

0 

L2 

0 

0 

0 

0  _ 

LL3- 

=  hCLQ  +  i^-.c  +  l2«s) 


(6.17) 


where  we  have  used  the  abbreviations   cos2<f>  =  C   and   sin2<f>  =  S.   As  with  initially 
unpolarized  radiation,  the  transmitted  radiation  is  completely  linearly  polarized  at  the 
orientation   <j)   of  the  polarizer  axis.   Equation  (6.17)  takes  on  an  especially  simple  form 
if  the  initial  ray  is  completely  linearly  polarized.   In  this  case,  as  we  have  seen,   |_ 
has  the  form 


L-Lr 


1 
cos2<|> 
sin2<J; 

0 


where  4>     is  the  polarization  direction.   Substituting  from  this  expression  L  =  L  'cos2^ 
and  L  =  L  «sin2^   in  eq.  (6.17)  we  find 
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|_'  =  ^L0"[l  +  cos2(<HC)] 


LQ«cos2((t.-^) 


1 

cos2<j> 

sin2<J> 

0 


(6.18) 


When  <$>   =   ty      the  polarizer  has  no  effect  on  the  ray  but  when   <J>  =  \p   ±  90°   the  ray  is  com- 
pletely extinguished.   Equation  (6.18)  is  known  as  Malus '  Law  and  expresses  the  familiar 
behavior  of  a  pair  of  ideal  polarizers. 

Before  going  on  to  consider  other  polarizing  components  we  should  remark  that  the 
discussion  above  can  be  turned  around  and  used  to  derive  the  Mueller  transmittance  matrix 
of  an  ideal  polarizer  [eq.  (6.15)].   The  argument  runs  as  follows:   We  define  an  ideal 
linear  polarizer  as  an  object  (1)  which  cannot  be  used  to  distinguish  circularly  polarized 
light  from  unpolarized  light,  (2)  whose  output  is  always  completely  linearly  polarized  in 
a  direction  parallel  to  its  polarization  axis,  and  (3)  which  attenuates  completely  linearly 
polarized  light  in  accordance  with  Malus'  Law  [eq.  (6.18)].   If  we  now  examine  eqs .  (6.7) 
we  see  that,  because  the  first  property  says  that  the  L  '   components  on  the  left  do  not 
depend  upon  the  value  of   L  ,   the  last  column  of  the  matrix  must  consist  entirely  of 
zeros.   The  second  property  says  that  the  output  of  the  polarizer  is  proportional  to 


1 
cos2<j 
sin2<j 

0 


where   <Ji   is  the  orientation  of  the  polarizer  axis.   Combining  this  with  the  third  property 
we  see  that  if  the  input  ray  has  a  radiance  vector 


L=Lr 


1 

cos2a 

sin2a 

0 


the  output  ray  will  be  described  by 


L1    =  L0-cos2(<j)-a) 


1 
cos2<{> 
sin2<J> 

0 


Therefore  the  Mueller  matrix  P(<f>)   for  an  ideal  linear  polarizer  must  satisfy 
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L'  =  *Q 


1  +  cos2<f>*cos2a  +  sin2<j>'sin2a 

cos2<j)  +  cos22<j>«cos2a  +  sin2<j>*cos2<f)'sin2a 

sin2cf)  +  sin2<}>*cos2<J>*cos2a  +  sin22<j>*sin2a 

0 


=  P(+)'Lr 


1 
cos2a 
sin2a 

0 


where  we  have  substituted   cos2(cf>-a)  =  %(1  +  cos2())-cos2a  +  sin2<}>'sin2a)  .   Comparing  this 
equation  with  eqs.  (6.7)  and  remembering  that  the .  last  column  of   P(4>)   must  be  zero  we 
immediately  obtain  eq.  (6.15).   This  derivation  illustrates  the  possibility  of  deducing 
the  form  of  a  Mueller  matrix  from  the  observed  experimental  behavior  of  the  object  it  is 
meant  to  describe.   We  will  later  use  this  approach  again  with  a  Mueller  matrix  form  which 
more  accurately  describes  real  (non-ideal)  polarizers. 

THE  IDEAL  LINEAR  RETARDER 

A  linear  retarder  is  an  optical  device  frequently  made  of  a  thin  plate  of  an  aniso- 
tropic crystal  such  as  calcite  or  quartz  with  its  faces  cut  parallel  to  the  crystal  optic 
axis.   Such  plates  have  the  (birefringence)  property  that  light  linearly  polarized  in  one 
direction,  called  the  fast  axis  of  the  retarder,  travels  with  a  higher  velocity  than  light 
linearly  polarized  in  a  direction  at  right  angles  to  this,  called  the  slow  axis.   A  mixture 
of  both  polarization  forms  entering  a  retarder  then  will  emerge  with  the  components  phase 
shifted  with  respect  to  one  another  and  consequently  with  its  polarization  form  changed. 
The  amount  of  phase  shift  is  called  the  retardance,   6.   The  retardance  is  related  to  the 

thickness,   a,   of  the  plate  and  the  indices  of  refraction  n^   and  n   for  the  fast  and 

f        s 

slow  polarizations,  respectively,  according  to 


2-Tra 


(nf  -  ns) 


[rad] 


(6.19) 


The  Mueller  transmittance  matrix  for  an  ideal  retarder  is: 


R(6,<f>)    = 


10  0  0 

0        C2  +  S2-cos6        S-C-(l  -  cos6)      -S-sin6 


0     S-C-  (1   -  cos6) 
0  S«sin6 


S2  +  C2-cos6 
-Osin6 


Osind 

cos6 


(6.20) 


where   C  =  cos2<J>,   S  =  sin2<f>   and   c|>   is  the  angle  between  the  fast  axis  of  the  retarder 
and  the  horizontal. 

Retarders  in  which  the  retardance  is  one-quarter  cycle   (6  =  tt/2)   or  one-half  cycle 
(5  =  it)   are  called  quarter-  and  half-wave  plates.   Quarter-wave  plates  are  useful  for 
converting  linearly  polarized  light  to  circularly  polarized  light,  and  vice-versa;  half- 
wave  plates,  for  rotating  the  direction  of  polarization  of  linearly  polarized  light.   If 
we  set   6  =  tt/2  in  eq.  (6.20)  we  obtain  the  Mueller  matrix  for  an  ideal  quarter-wave  plate: 
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1 

0 

0 

0 

0 

c2 

s-c 

_c 

0 

s-c 

s2 

c 

0 

s 

-c 

0 

Q($)  =  R(tt/2,4>)  =   „  ...   o2   „  (6.21) 


DETERMINING  STOKES  COMPONENTS  with  IDEAL  INSTRUMENTS 

With  an  ideal  linear  polarizer  and  an  ideal  quarter-wave  plate  we  are  now  in  a  posi- 
tion to  measure,  in  principle,  the  four  Stokes  components  of  an  arbitrary  light  ray.   We 
pass  the  ray  through  the  quarter-wave  plate  and  then  through  the  polarizer  and  measure  the 
output  with  a  polarization-indifferent  detector.   The  Mueller  transmittance  matrix  for  the 
combination  of  quarter-wave  plate  at  angle   <f>  =  a   and  polarizer  at  angle   <j>  =  3   is 

A  =  P(3)-Q(a) 


1  cos23 

cos23  cos223 

sin23  sin23-cos2f 

0  0 


sin23  0 

sin23-cos23  0 

sin223  0 

0  0 


10  0  0 

0    cos  2a  sin2occos2a  -sin2a 

0   sin2a*cos2a  sin  2a  cos2a 

0     sin2a  -cos2a  0 


1       cos2a#cos2(3-a)        sin2a*cos2(3-a)  sin2(3-a) 

cos23  cos23*cos2a'Cos2(3-a)  cos23#sin2a*cos2 (3-a)  cos23#sin2 (3-a) 

sin23  sin23*cos2a'cos2(3-a)  sin23*sin2a-cos2(3-a)  sin23*sin2(3-a) 

0                                0                                                 0  0 


(6.22)1 


If   a  ray  with  spectral-radiance  vector 


LL3J 


passes  through  this  combination,  the  spectral  radiance  vector   [_'   of  the  output  ray  will 
be  given  by 

1 
cos23 
sin26 

0 


L'  =  A*L  =  ^[L0  +  L  -cos2a-cos2(3-a)  +  L2'sin2a-cos2(3-a)  +  L3"sin2  (3-a)] 


(6.23) 


Obviously,  unless  one  can  find  simplifying  trigonometric  identities  a  long  train  of 
polarizing  optics  will  quickly  lead  to  an  almost  unmanageable  problem.   The  operator  form 
of  Mueller  matrix  manipulation  described  by  Priebe  [9]  permits  many  simplifications  to  be 
performed  before  the  matrix  elements  become  too  unwieldy.   His  technique  is  especially 
useful  for  weakly  (linearly)  polarizing  components,  small  retardances ,  and  small  angles 
of  rotation. 
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Since  we  are  only  interested  in  the  response  of  a  polarization- indifferent  detector  to  this 
ray  we  actually  need  only  the  L  '   component  of   [_'   which  is 

LQ'  =  %[LQ  +  L1-cos2cccos2(6-a)  +  L  • sin2a-cos2(B-a)  +  L  -sin2(3-a)]  .     (6.24) 

Had  we  been  interested  in  passing  the  output  into  or  through  any  further  polarization- 
sensitive  components  we  would,  in  general,  have  needed  all  four  components  of   |_' .   From 
eq.  (6.24)  we  see  that  if  we  make  four  measurements  of   L  '   at  four  combinations  of   a 
and   3,   we  can  solve  the  four  resulting  equations  for  the  Stokes  components   L_,   L  , 
L„,   and   L   of  the  original  ray.   A  particularly  simple  solution  is  obtained  when  a  and 
B   are  set  at   0°,  45°,  and  90°   as  follows: 

LJ  (a  =  0,  6  =  0)  =  lg(L0  +  L1) 

L'  (a  =  0,  6  =  90°)  =  H(L     -  L. ) 

(6.25) 
LQ'  (a  =  45°,  B  =  45°)  =  %(LQ  +  L2) 

Lj  (a  =  0,  B  =  45°)  =  %(LQ  +  1^). 

From  these  four  measured  values  of   L  '   the  polarization  components  of  the  original  ray 
can  be  obtained  at  once. 

In  general  the  result  of  passing  a  ray  through  an  ideal  retarder  followed  by  an  ideal 
linear  polarizer  is  a  ray  whose  spectral  radiance  is  a  sum  of  the  Stokes  components  of  the 
original  ray  each  multiplied  by  sines  and  cosines  of  the  orientation  angles  and  the 
retardance.1  Such  a  relationship  readily  lends  itself  to  Fourier  analysis,  and  many  auto- 
matic techniques  for  measuring  some  or  all  of  the  Stokes  components  are  based  on  this  idea 
[10,11]. 

It  is  possible  to  deduce  relative  values  of   L  ,   L  ,   and   L   by  measuring  the  pair 
of  angular  settings   a   and   B   of  the  quarter-wave  plate  and  polarizer  which  simultane- 
ously minimize  L  '  .   This  is  a  common  technique  in  ellipsometry,  for  example,  because  sl 
a  determination  requires  only  monotonicity  of  detector  response  —  not  linearity  or,  indeed, 
not  even  significe 
adjusting   a   and 


ously  minimize  L  '  .   This  is  a  common  technique  in  ellipsometry,  for  example,  because  such 

re 
not  even  significant  stability.   The  minimum  value  of   L  '   which  can  be  achieved  by 


-J  -  ^0  -  <h2  +  S2  +  h^\ 


This  occurs  when 


tan2a  =  L„/L 
and  (6.26) 

tan2(B-a)  =  1^/(1^  +  L^)2, 


1The  complete  expression  is  given  in  eq.  (6.41) 
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as  can  be  obtained  by  differentiation  of  eq.  (6.24).   From  a  measurement  of  the  angles   a 
and   B-a   at  this  minimum1  then  we  can  calculate  the  relative  values  of  the  polarization 
components  from  eq.  (6.26): 


L, /L  =  cos2a-cos2  (6-ot) 

1  P 

L„/L  =  sin2a'cos2(3-a) 

2  P 

L_/L   =  sin2(g-a), 

3  P 


(6.27) 


where  L  =  (L,2  +  L2  +  L2)   is  the  polarized  component  of   L. .   Unfortunately  the 
p     1     2     3  0 

determination  of   L.   and   L   cannot  be  carried  out  by  measurements  of  angles  only  but 
Op 

requires  flux  measurements  with  the  ubiquitous  problems  of  temporal,  positional,  direc- 
tional, and  wavelength  responsivity  and  of  detector  linearity  which  characterize  and 
bedevil  radiometric  measurements. 

An  alternative  way  to  measure  relative  values  of   L  ,   L  ,   and   L   is  by  using  a 
variable  retarder  —  usually  called  a  compensator .   This  is  a  retarder  composed  essentially 
of  two  wedge-shaped  retarders  arranged  so  that  sliding  them  across  one  another  changes  the 
total  thickness  of  the  composite  block  and  hence  the  retardance  of  the  device.   By  simul- 
taneously adjusting  the  retardance  and  angular  orientation  of  a  perfect  compensator, 
arbitrarily  polarized  light  can  be  converted  to  any  other  arbitrary  polarization  state 
which  possesses  the  same  degree  of  polarization.   Thus,  we  may  pass  radiation,  whose  Stokes 
components  are  to  be  determined,  through  a  compensator  and  then  through  a  fixed  perfect 
polarizer  and  adjust  the  compensator  so  that  the  polarizer  is  able  to  extinguish  the  polar- 
ized component  of  the  light.   From  the  retardance  and  orientation  settings  of  the  compen- 
sator when  this  adjustment  has  been  achieved  the  Stokes  component  values  can  be  computed. 
Taking  for  simplicity  a  polarizer  with  its  polarization  axis  horizontal  we  have  from 
eqs.  (6.15)  and  (6.20) 


P(0°).R(6,4>).|_  =  h 


1  1  o  o" 

110   0 

0   0  0  0 

0   0   0   0_ 

10  0  0 

0      C2  +  S2-cos6  S«0(l-cos6)  -S-sin^ 

0      S-0(l-cos5)  S2  +  C2*cos6  C-sin<! 

0             S«sin6  -Osin6  cos6 


J0 


.S-. 


The  output  radiance  is 

L»  =  %  [l  +  L1-(C2  +  S2-cos5)  +  L2-S-C-(l-cos6)  -  I^-S-sinfi], 


*In  the  literature  of  ellipsometry  the  angles  a  and  x  =  &-<*  defined  by  eq.  (6.26)  are 
known  respectively  as  the  azimuth  and  ellipticity  of  the  polarized  component.  See  refer- 
ence [6]  for  a  review  of  ellipsometry. 
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By  differentiation  we  find  that   L  '   has  a  minimum  value  of 


LQ'  -h   [Lq-CL,2  +  1^  +L/)1*]  =>2(L0-Lp) 


at  values  of  the  orientation  angle   <f>   and  retardance   6   which  satisfy: 

Ln/L  =  -(cos22<j>  +  sin22(j)'cos6) 

1  P 

L_/L     =   -sin2d>>cos2<j>- (1   -cos6) 

2  P 

L_/L  =  sin2d)'sin6. 

3  P 


(6.28) 


These  ellipsometric  techniques  can  be  extremely  accurate  and  are  used  extensively  in 
studying  surfaces  and  surface  films  by  the  analysis  of  the  relative  polarization  form  of 
reflected  light  [5,6].   Note  that  these  techniques  require  an  accurate  means  of  sensing  a 
minimum  flux  level  and  that   L   and  the  more  important  component   L„   remain  undetermined, 

THE  IDEAL  CIRCULAR  POLARIZER 

In  eqs.  (6.25)  we  saw  how  a  quarter-wave  plate  followed  by  an  ideal  linear  polarizer 
with  its  polarization  axis  at  45°  to  the  fast  axis  of  the  retarder  could  be  used  to  gather 
information  about  the  Stokes  component   L  ,   which  is  the  circular  polarization  component. 
In  this  case  the  Mueller  transmittance  matrix  [eq.  (6.22)]  reduces  to: 


Acir  =  P(3)'Q(a  =  6+45°)  = 


1  0  0  ±1 

C  0  0  ±c 

S  0  0  ±s 

0  0  0  0 


(6.29) 


where   C  =  cos  23,   S  =  sin  23,   3   describes  the  orientation  of  the  linear  polarizer  and 
the  +  or   -   signs  are  chosen  according  to  the  sign  of   6  —  a  =  +45°.   Arbitrarily 
polarized  light  passing  through  this  combination  will  emerge  with  Stokes  vector 


L'  =  A 


cir 


~v 

1 

Ll 

L2 

"  hiL0   +-  V  ' 

cos2 
sin2 

_L3_ 

0 

(6.30) 


Suppose,  for  example,  that  the  incoming  light  is  completely  circularly  polarized  so  that 

L_  =  ±L„   and  that  its  handedness  is  opposite  to  that  of  A  .  .   We  will  then  observe 
0     3  cir 

[_'  =  0,   i.e.  complete  extinction.   Notice  that  although  a  measurement  of  the  spectral 
radiance  of  the  emergent  ray  yields  information  primarily  about  the  circular  polarization 
component   L„   the  emergent  ray  itself  is  linearly  polarized.   In  order  to  generate  cir- 
cularly polarized  light  we  must  reverse  the  retarder-polarizer  combination:   the  light  must 
first  pass  through  the  linear  polarizer  and  then  through  the  quarter-wave  plate.   Taken  in 
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this  order  the  combination  has  the  Mueller  matrix 


Pcir  =  Q(3  ±45°).P(3)  =  h 


1  C  S  0 

0  0  0  0 

0  0  0  0 

±1  +C  ±S  0 


(6.31) 


where   C  =  cos2g,   S  =  sin23   and   6   is  the  orientation  of  the  linear  polarizer.   The 
quarter-wave  plate  is  at  the  angle  a  =  8  ±45°.   Now  if  arbitrarily  polarized  light  passes 
through  this  combination  the  result  is 


L'  -P. 


cir 


LL3 


kCLQ   +  L  cos28  +  L  sin2g) 


1 

0 

0 

±1 


(6.32) 


which  represents  pure  circularly  polarized  light.   In  addition  to  illustrating  the  mathe- 
matical quirk  that  the  product  of  two  matrices  may  depend  upon  the  order  in  which  they  are 

taken,  the  difference  between  A  .    and   P  .    [(or  eqs.  (6.29)  and  (6.31)1  is  important 

cir        cir   ltv/  t- 

experimentally  since  commercial  circular  polarizers  are  usually  made  by  cementing  a 
quarter-wave  plate  to  a  linear  polarizer  with  their  axes  45°  apart.   Such  a  device  will 
act  as  a  circular  polarizer  only  if  light  enters  from  the  linear  polarizer  side  and  will 
act  as  an  analyzer  of  circularly  polarized  light  only  if  the  light  enters  from  the  quarter- 
wave  plate  side. 

An  ideal  circular  polarizer  without  these  idiosyncracies  can  be  made  by  adding  a 
second  quarter -wave  plate  to  make  a  sandwich  with  the  linear  polarizer  in  the  middle. 
Each  element  must  be  turned  45°  from  the  previous  one  to  give 


cir 


=  Qfcj,  ±90°)-P(^  ±45°)-QO) 


(6.33) 


where  Q(<j>)   is  the  Mueller  matrix  [eq.  (6.21)]  for  a  quarter -wave  plate  and   P(t)))   is  that 
for  a  linear  polarizer  [eq.  (6.15)].   If  we  multiply  out  the  matrices  we  find 


P  .  =  h 

1  cir 


1   0  0  ±1 

0   0  0  0 

0   0  0  0 

±10  0  1 


(6.34) 


Plus  signs  in  eq.  (6.33)  lead  to  plus  signs  in  eq.  (6.34)  and  a  right-circular  polarizer. 
Minus  signs  throughout  give  a  left-circular  polarizer.   This  circular  polarizer  now  has 
properties  analogous  to  those  we  associate  with  the  ideal  linear  polarizer:   it  can  be 
used  as  polarizer  or  analyzer  from  either  side;  it  will  pass  circularly  polarized  light 
of  the  proper  handedness  but  will  extinguish  any  component  of  the  opposite  handedness;  and 
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it  cannot  distinguish  between  linearly  polarized  light  and  unpolarized  light. 

RADIOMETER  CALIBRATION 

Equation  (6.32)  is  a  special  case  of 

[_'  =  Q(a)-P(0)-L 

1  cos2g  sin23  0 

cos2a*cos2 (a-g)  cos2a-cos2g«cos2(a-g)  cos2a*sin2g'cos2(a-g)  0 

sin2a*cos2(a-g)  sin2a«cos2g-cos2(a-g)  sin2a-sin2g«cos2(a-g)  0 

sin2(a-g)        cos2g  sin2(a-g)        sin2g-sin2(a-g)  0 


(6.35) 


%(L  +  L  -cos2g  +  L  «sin2g) 


cos2occos2 (a-g) 

sin2a*cos2 (a-g) 

sin2(a-g) 


which  can  be  obtained  by  multiplying  the  matrices  of  eqs.  (6.21)  and  (6.15).   This 
describes  the  behavior  of  a  polarizer  at  orientation  angle   g   followed  by  a  quarter-wave 
plate  at  orientation  angle  a.   By  the  proper  orientation  of  the  quarter-wave  plate  and 
polarizer  any  arbitrary  state  of  polarization  of  the  emerging  beam  can  be  achieved.   If 
g   =   a,   for  example,  the  output  is  linearly  polarized  at  the  orientation  a,   and  if 
g  =  a  ±45°  the  output  is  right  (+)  or  left  (-)  circularly  polarized.   Such  a  beam  can  be 
used  to  measure  Mueller  transmittance  matrix  elements  and  to  calibrate  a  radiometer. 

The  complete  Mueller  transmittance  matrix  for  an  optical  component  must  be  known  if 
the  polarization  state  of  a  ray  after  passing  through  this  component  is  wanted.   However, 
in  many  cases  all  that  is  needed  is  a  prediction  of  the  response  of  a  complete  instrument 
to  the  input  radiation.   This  is  a  simpler  problem  because  the  instrument  response  depends 
only  upon  the  spectral-radiance  component   L  '   which  triggers  the  final  electronic, 
thermal,  or  chemical  response.1   Here,  only  the  top  row  of  the  Mueller  matrix  of  the 
instrument  is  required.   That  is,  only  the  first  of  eqs.  (6.7)  is  needed.   We  can  write 

dS  =  (R0().L0  +  R01.Ll  +  RQ2.L2  +  R03.L3)  -dG-dA 

where  dS   is  the  element  of  radiometer  signal  output  associated  with  the  ray  whose  Stokes 
spectral-radiance  components  are  L  ,   L  ,   L  ,   and   L   and  whose  elements  of  throughput 


■^Some  detectors  may  respond  to  properties  of  the  electromagnetic  radiation  other  than 
power  —  angular  momentum,  for  example.   It  appears  that  such  detectors  can  still  be 
treated  in  a  purely  formal  way  as  we  describe  here,  although  the  resulting  responsivities, 
R   ,   may  not  be  proportional  to  the  Mueller  matrix  elements  of  any  physically  realizable 
optical  component. 
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and  wavelength  are  dO  =  dwcosO'dA  and  dA,  respectively.   The  coefficients   Rnn>   Rni  ' 
etc.,  are  the  first-row  Mueller  transmittance-matrix  elements  of  the  optical  part  of  the 
radiometer,  including  the  optical  polarization  characteristics  of  the  detector-  or 
transducer-element  responsivity,  multiplied  by  the  polarization-indifferent  responsivity 
of  that  element  in  combination  with  any  signal-processing  and  display  components  of  the 
instrument.    The  total  response  of  the  instrument  is  the  integral  of  this  expression  over 
all  rays  which  can  reach  the  photo-sensitive  surface. 


S  = 


AA 


AA 


(RQ0-L0  +  RQ1'L  +  Rq2'L2  +  R03'L  )-dwcose«dA"dA  (6.36) 

Acj 


where   AA,   AA,   and   Aw   are  the  acceptance  intervals  of  the  radiometer.   If  these  accept- 
ance intervals  are  small  enough  that  the  Stokes  radiance  vector  can  be  considered  uniform 
over  all  of  them  then2 

S  =  *00'L0  +  *0l'Ll  +  *02-L2  +  *03'L3  (6-36a) 

where  Rn .    = 
Oi 


Rn_.  'du'cosO  *dA'dA.   These  four  responsivity  factors  RnA      can  be 

AA;  Aw 


AAJ 

determined  by  observing  the  response  of  the  instrument  to  four  sufficiently  different, 
known  polarization  states  of  a  test  beam  with  fixed  geometry  and  spectral  distribution. 
We  thus  obtain  four  equations  of  the  form  of  eq.  (6.36a)  to  be  solved  for  the  four  respon- 
sivity factors  Rn-'      The  L.   are  the  Stokes  components  of  the  input-beam  radiance  pro- 
duced, for  example,  as  described  above,  following  eq.  (6.35),  and  measured,  if  necessary, 
as  described  in  connection  with  eq.  (6.24)  using  a  polarization-indifferent  detector  and 
a  second  quarter-wave  plate  and  linear  polarizer. 


If  the  acceptance  intervals  of  the  instrument  are  too  large  to  assume  uniform  radiance, 
then  it  may  be  that  the  coefficients   R    in  eq.  (6.36)  can  be  considered  constant,  inde- 
pendent of  wavelength  and  of  ray  position  and  direction.   In  this  case 


Must  as  the  electrical  engineer  replaces  an  actual  generator  by  an  equivalent  circuit 
consisting  of  either  an  ideal  voltage  generator  and  series  impedance  or  an  ideal  current 
generator  and  shunt  impedance,  it  is  convenient,  analytically,  to  replace  the  actual 
detector  or  transducer  element  by  an  equivalent  combination  of  (1)  a  pure  transmittance 
(whose  Mueller  transmittance  matrix  can  be  combined  with  that  of  the  other  optical  ele- 
ments by  matrix  multiplication)  followed  by  (2)  a  polarization-indifferent  detector  ele- 
ment that,  in  combination  with  the  following  signal-processing  and  display  components, 
determines  the  magnitude  of  the  overall  responsivity. 

As  in  earlier  chapters  we  assume  linearity  throughout  this  discussion.   That  is,   R_ . 
may  depend  upon  ray  position,  direction  and  wavelength  but  not  upon  Stokes  spectral- 
radiance  components. 
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where   $ . 

1 


R   •$   +  R   •$   +  R   •$   +  R   •$  (f,    16M 

00  0    01  1    02  2    03  3  ko.JOD; 


L . -dwcose  *dA-dA.   Again  the  four  responsivity  coefficients  can  be 

'  L\'  AAJ  Aw 


determined  from  observations  of  the  output  response   S   to  four  beams  with  known,  different, 
Stokes  spectral-radiance-vector  distributions.   A  common  way  of  achieving  approximately 
isotropic  (directionally  uniform)  responsivity,  at  least  over  the  acceptance  solid  angle 
Aw,   is,  as  we  have  seen  in  Chapter  5  [12],  to  introduce  a  diffuser  as  the  first  component 
of  the  radiometer.   In  this  case,  since  such  a  diffuser  is  also  a  depolarizer,  we  have  the 
further  simplification:   R  ..  =  R   =  R   =  0.   The  quantity   <3>   is  the  total  radiant  flux 
in  the  wavelength  pass  band  of  the  radiometer. 

If  neither  the  R„ .   nor  the  L.   are  uniform, the  dependence  of  the   R„ .   upon  ray 
Oi  l  Oi 

direction,  position  and  wavelength  can  still  be  determined  by  probing  the  instrument  with 
well-defined,  narrow  pencils  of  rays  with  uniform  Stokes  spectral-radiance  vector.  The 
general  usefulness  of  such  a  radiometer,  however,  is  questionable  since  its  response  is  a 
sum  of  weighted  integrals  of  the  spectral  radiance  distribution  falling  upon  it.  For  com- 
parison of  sources  with  identical  relative  spectral  radiance  distribution  over  direction, 
position  and  wavelength,  of  course,  such  a  radiometer  might  be  satisfactory  and  its  cali- 
bration in  this  case  again  calls  for  only  four  measurements. 

THE  ROTATION  MATRICES 

If  we  look  back  at  eqs.  (6.15)  and  (6.21)  we  notice  a  similarity  in  the  way  the  angle 
<j)   enters  into  the  matrices  for  a  polarizer  and  for  a  quarter-wave  plate.   This  is  not 
accidental.   It  arises  in  the  following  way:   suppose  we  somehow  were  given  numerical 
Mueller  propagance-matrix-element  values  for  an  optical  component  at  a  particular  orienta- 
tion, say   <j)  =  0.   How  can  we  generalize  this  to  other  orientations,   4>  4-   0?   We  do  it  in 
a  kind  of  backward  way  by  first  finding  the  Stokes  components  of  an  entering  ray  in  a 
coordinate  frame  which  has  been  rotated  through  the  angle   cj>   so  as  to  match  the  orienta- 
tion for  which  we  know  the  Mueller  matrix  of  the  (rotated)  optical  component.   Then  we  can 
multiply  this  Stokes  vector  by  the  Mueller  matrix  and  finally  rotate  the  coordinate  frame 
back  to  the  original  "horizontal"  orientation. 

Consider  an  arbitrarily  polarized  ray  described  by 

I 


J0 

Ll 

L 


2 
LL3 

If  we  rotate  our  coordinate  frame,  only  the  linear  polarization  components,   L   and  L  , 
will  be  affected  because  these  are  the  only  components  whose  definitions  require  specifica- 
tion of  a  reference,  horizontal  direction.   So,  thinking  of   L   and  L   for  the  moment  as 
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orthogonal  components  of  a  two-component  vector  (see  fig.  6.2)  we  can  easily  find  the  new 

components  of  this  vector  in  a  rotated  frame  of  reference.   In  real  space  we  know  that  if 

we  rotate  the  "horizontal"  direction  by  45°  we  will  have  rotated  the  direction  of  L   into 
that  of   ±L  .   In  our  two-dimensional  vector  space,  however,  the  angle  between  L   and   L 

is  90°;  therefore  a  rotation  through  the  angle   <j>   in  real  space  corresponds  to  a  rotation 
through  the  angle   2<J>   in  our   (L  ,L  )   vector  space.   We  see  from  figure  6.2  that  such  a 
rotation  of  our  coordinate  axes  will  leave  us  with  new  components 

L'  =   L  •cos2<()  +  L  'Sin2<}> 

and 

L'  =  -L  'Sin2cj)  +  L  •cos2<j>. 

This  result,  combined  with 

LJ  =L0 

and 

L3  =L3 
is  very  suggestive  of  eq.  (6.7):   that  is,  by  defining  a  matrix  M(<j>)   we  can  write 


L"  =  M(*)-L- 


10  0  0 

0        cos2<(>  sin2<J)  0 

0      -sin2<}>  cos2cf>  0 

0  0  0  1 


(6.37) 


to  describe  this  rotation.   The  procedure  for  handling  an  optical  component  turned  at  an 
angle   (f>   is  now  fairly  clear.   Rotating  the  coordinates  of     through  the  angle   <j>   so 
that  "horizontal"  agrees  with  the  optical-component  orientation,  applying  its  Mueller 
propagance  matrix,  and  rotating  again  through   -<j>   to  get  back  to  the  real  horizontal  is 
equivalent  to  the  following 


L"  = 


10  0  0 

0       cos(-2<|>)  sin(-2(j))  0 

0      -sin(-2tj))  cos(-2<fi)  0 

0  0  0  1 


10  0  0 

0        cos2<j>  sin2<j>  0 

0      -sin2cj>  cos2<(>  0 

0  0  0  1 


(6.38) 


in  which   ;   is  the  Mueller  matrix.   The  matrix  M(<}>)   defined  in  eq.  (6.37)  is  known  as 
a  rotation  matrix.   If  the  multiplication  of  the  three  4x4  matrices  in  eq.  (6.38)  is 
carried  out  we  obtain  the  equivalent  Mueller  matrix  evaluated  at  the  angle  <j>.   That  is, 
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L  '  =  OD  =  OA  +  AD  =  OA  +  BC 


=  L  'cos2(j)  +  L  «sin2q 


L2'  =  DP  =  CP  -  DC  =  CP  -  AB 


=  L  'COs2<f>  -  L  «sin2c 


Figure  6.2.    The  components  of  £   in  the  1^,  L2   and   L  *   L  '   coordinate  systems. 
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T«>>  =  M(-*)-T(o)-M<*) 


oo 


T10C   T20S 


T10*S  +  T20-C 


"30 


T0rC  "  T02'S        T0l"S  +  T02*C 


Tn-A 


T21  +  B 


12 


T22+A 


T3l'C   "   T32*S        T3l'S   +  T32*C 


"03 


T13*C  "'  T23*S 


T13.S  +  T23.C 
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(6.39) 


where   A  =  (T   -  T^-S2  +  (T  2  +  T^-S-C 

B  =  (j        -  T22)-S-C  -  (T   +  T21)-S2 

S   =   sin2cj> 
C  =   cos2<f> 

and  Tnn»   Tm  '   Tn?'   etc«>  are  taken  to  be  the  matrix  elements  [eq.  (6.10)]  at   <}>  =  0. 
Thus,  in  their  dependence  upon  the  orientation  angle  $,      all  Mueller  matrices  must 
resemble  the  matrix  product  of  eq.  (6.39).   This  multiplication  has  already  been  performed 
for  most  of  the  matrices  we  discuss  in  this  chapter  and  accounts  for  the  angular  orienta- 
tion dependence  built  into  them.   Many  Mueller  matrices  in  the  literature,  however,  are 
written  only  for  a  special  angle,  usually   <j>  =  0,   and  require  the  transformation  of 
eq.  (6.39)  before  they  can  be  used  at  other  orientations. 

In  addition  to  being  useful  in  providing  the  dependence  of  the  transmittance  of 
optical  components  upon  orientation,  the  rotation  matrix  M(<j))   of  eq.  (6.37)  is  the 
Mueller  transmittance  matrix  for  an  ideal  circular  retarder.   Circular  retarders  are  made 
of  materials  which  exhibit  the  property  called  optical  activity  or  circular  birefringence. 
Crystal  quartz  and  solutions  of  dextrose  are  examples  of  such  materials.   Linearly  polar- 
ized light  entering  a  circular  retarder  will  emerge  with  its  direction  of  polarization 
rotated  by  the  angle   $   appearing  in  ^\(<$>)  .   This  angle  is  the  circular  retardance  and  is 
related  to  the  wavelength  X     and  the  thickness   a  of  the  material  by 

<p   =  -~-   (nr  -  n£)   [rad] 


[compare  with  eq.  (6.19)]  where  n   is  the  refractive  index  for  right-handed  circularly 
polarized  light  at  wavelength   A   and   n   is  for  left-handed.   For  solutions  of  an  opti- 
cally active  material  n  -  n.   will  be  proportional  to  the  concentration  of  solute  so  that 
measurements  of  the  angular  rotation  of  the  direction  of  polarization  for  a  fixed  path 
length  provide  a  convenient  means  of  measuring  solute  concentrations.   This  technique  has 
been  developed  to  a  high  degree  of  accuracy  in  instruments  called  saccharimeters  or 
polarimeters. 
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A  different  kind  of  rotation  problem  arises  if  we  have  the  Mueller  matrix  elements  of 
an  optical  component  or  train  of  components  and  ask  what  happens  if  the  radiation  passes 
in  the  opposite  direction  through  the  system.   How  does  the  Mueller  matrix  change?   The 
answer  is  that  for  an  original  matrix  of  elements   T..   given  by  eq.  (6.10)  the  matrix  for 
the  reversed  optical  system  is 


■(r) 


00 


01 


02 


03 


'10 


-T 


20 


'30 


lll 


-T 


21 


'31 


-T 


12 


22 


-T 


"13 


-T 


23 


32 


■33 


(6.40) 


the  rows  and  columns  have  been  interchanged1  and  the  signs  of  the  matrix  elements  are 
changed  once  for  every  subscript  2  which  appears.   The  interchange  of  rows  and  columns  is 
understandable  because  the  individual  columns  operate  on  the  individual  components  of  the 
Stokes  vector  of  the  input  radiation  while  the  rows  generate  the  components  of  the  Stokes 
vector  of  the  output  radiation.   So  if  we  interchange  the  roles  of  input  and  output  it 
isn't  surprising  that  the  roles  of  the  rows  and  columns  must  be  interchanged  too.   The 
change  of  sign  associated  with  the  subscript  2  is  caused  by  the  transformation  of  the  +45° 
polarization  component  into  its  complementary  +135°  component  when  an  object  is  reversed.2 
Imagine  a  plate  containing  milled  slots  in  horizontal,  vertical,  and  45°  directions  and 
also  a  hole  tapped  for  a  right-handed  screw.   Now  if  the  plate  is  viewed  from  the  back  the 
horizontal  and  vertical  slots  and  the  tapped  hole  will  appear  unchanged  but  the  45°  slot 
will  look  like  a  135°  slot.   Thus,  the  +45°  and  +135°  directions  change  roles  upon  rever- 
sal, and  this  requires  the  sign  changes  shown  in  eq.  (6.40). 

REAL  OPTICAL  COMPONENTS 

So  far  in  this  chapter  we  have  seen  that,  if  we  have  a  source  of  unpolarized  light,  a 
polarization-indifferent  detector,  a  perfect  linear  polarizer,  and  a  perfect  quarter-wave 
plate,  we  can  determine  the  state  of  polarization  of  any  arbitrary  ray  and  by  varying  the 
polarization  of  such  a  ray  we  can  use  it  to  characterize  or  calibrate  the  response  of  an 
instrument  to  polarized  light.   Unfortunately  real  polarizers  and  quarter-wave  plates 
exhibit  Mueller  transmittance  matrices  which  differ  somewhat  from  those  given  in 


*If  the  rows  and  columns  of  a  matrix  are  interchanged  the  result  is  known  as  the  transpose 
of  the  original  matrix. 

2It  is  not  obvious  that  the  relationship  between  the  matrix  for  the  combination  of  a  polar- 
izer followed  by  a  quarter -wave  plate  in  eq.  (6.35)  and  its  reverse  in  eq.  (6.22)  satisfies 
eq.  (6.40)   This  is  because  the  angles  used  in  eq.  (6.35)  are  the  supplements  of  those  in 
eq.  (6.22);  when  an  optical  component  at  orientation   <f>   is  reversed  it  will  find  itself  at 
orientation  (180°  -  <j>)  .   If  the  matrices  in  eq.  (6.35)  and  (6.22)  are  both  expressed  in 
terms  of  the  same  angles  it  will  be  found  that  eq.  (6.40)  is  obeyed. 
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eqs.  (6.15)  and  (6.21).   High  quality  crystal  polarizers  and  retarders  will  probably 
approach  the  matrices  of  eqs.  (6.15)  and  (6.20)  closely  except  for  the  effects  of  surface 
reflections.   So,  for  many  radiometric  applications,  where  polarization  complications  are 
no  more  than  small  corrections  to  measurements  often  made  with  complete  neglect  of  polari- 
zation, such  optical  elements  can  be  used  and  the  Mueller  matrices  of  eqs.  (6.15)  and  (6.20) 
will  be  adequate  to  account  for  their  behavior  (to  within  10%  or  so) . 

Regardless  of  whether  a  polarizer  and  retarder  are  perfect,  if  their  Mueller  matrices 
are  known  the  Stokes  components  of  a  ray  can  still  generally  be  measured  essentially  as 
described  in  connection  with  eq.  (6.24).   Four  or  more  combinations  of  angular  settings  of 
the  polarizer  and  retarder  are  needed  as  before.   The  difference  is  that  the  Mueller  matrix 
of  the  combination  of  imperfect  elements  will  differ  from  eq.  (6.22)  so  that  eq.  (6.24) 
will  not  be  valid.   However,  the  analogue  to  eq .  (6.24)  will  still  be  linear  in  the  com- 
ponents L„,   L..  ,   L„,   and  L   and  the  four  measurements  at  suitable  angular  settings 
will  permit  these  components  to  be  evaluated.   In  the  simplest  case  the  retarder  will  not 
be  exactly  a  quarter-wave  plate  due,  for  example,  to  its  use  at  a  wavelength  other  than  the 
one  for  which  it  was  designed.   For  this  case  eq.  (6.24)  becomes  [using  eq.  (6.20)  in  (6.22) 
instead  of  (6.21)]: 


{ 


L  +  L  • [cos  2a*cos  2(B-a)  -  sin  2a*sin  2 (6-a) #cos6] 

+  L  -[sin  2a*cos  2(g-a)  +  cos  2a*sin  2(g-a)'cos6]        (6.41) 


+  L  • [sin  2(3-a)-sin6]| 


Obviously  if   6   is  known,  from  eq.  (6.19)  for  example,  this  equation  is  no  more  difficult 
to  work  with  than  eq.  (6.24). 

The  real  problem  arises  when  one  doesn't  know  the  Mueller  matrices  of  his  polarizers 
or  retarders  and  has  no  sources  of  polarized  light  with  known  Stokes  components  which  can 
be  used  to  calibrate  them.   One  can,  with  considerable  confidence,  produce  unpolarized 
light  and  a  polarization- indifferent  detector  (using  diffusers,  such  as  integrating 
spheres  [13],  for  instance).   With  such  a  source  and  detector,  one  can  determine  only 
element   T    of  the  Mueller  matrix  of  an  optical  component  or  combination  of  components. 
However,  as  we  have  seen  in  eq.  (6.41),  for  example,  when  polarizing  components  are  com- 
bined the  transmittance  to  unpolarized  radiation  of  the  combination  depends  upon  additional 
matrix  elements  of  the  individual  components.   Thus,  we  might  expect  that  measurements  of 
the  transmittances  of  all  combinations  of  several  polarizing  components  as  a  function  of 
their  orientation  angles  can  yield  considerable  information  about  the  matrix  elements  of 
the  individual  components.   Using  this  idea  we  show  in  Appendix  7  how  almost  any  Mueller 
matrix  can  be  measured  using  nothing  more  than  an  unpolarized  light  source  and  a  polariza- 
tion-indifferent detector.   In  the  following  paragraphs  we  show  how  the  same  technique  can 
be  applied  to  measure  a  simple  matrix  which  approximately  describes  the  common  Polaroid 
type  of  polarizer. 
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Our  procedure  simultaneously  measures  three  linear-polarizing  filters.1   We  assume 
that  each  is  mounted  in  such  a  way  that  it  can  be  reproducibly  inserted  and  removed  from 
the  unpolarized  measurement  radiation  beam  and  that  it  can  be  rotated  about  the  beam  axis, 
with  its  angular  orientation  read  from  a  scale  fixed  to  the  filter  mount.   Real  polarizers, 
especially  sheet  polarizers,  can  be  expected  to  exhibit  a  small  amount  of  retardance. 
Likewise,  real  retarders  may  be  slightly  linearly  polarizing.   Therefore,  for  the  Mueller 
transmittance  matrix  of  such  real  components  we  assume  a  form  which  is  general  enough  to 
include  both  functions: 


F(*)   = 


s  d«C  d-S  0 

d-C  s-C2  +  p-S2  (s-p)-S-C  -q-S 

d-S  (s-p)-S-C  s-S2  +  p-C2  q-C 

0  q-S  -q«C  p 


(6.42) 


In  this  expression 


the  filter  mount.   C  =  cos2(tj)-<t>  )   and   S  =  sin2(cf>-<f)  ),   where   cf> 


is  the  orientation  of  the  filter  as  read  from  the  scale  attached  to 

is  the  orientation  of 

the  polarization  axis  of  the  polarizing  filter.   The  angle   <J>~   and  the  parameters   s,  d, 
p,  and  q   are  to  be  determined  by  experiment.   The  form  of   FC^) >   °f  course,  includes  the 
<j> -dependence  of  eq.  (6.39). 

The  parameters   s,  d,  p,  and  q   represent  a  model  for  dichroic  polarizers  of  the 
Polaroid  type  with  allowance  for  some  retarder  behavior.2   In  this  model 


s  =  h(K    +   kj 


Vv 


k  *cosd 


where  the  retardance  of  the  filter  is   6 
principle  transmittances  of  the  filter.3 


s(k,  -  kj 


V 


k  *k  "sin6, 


and  the  quantities   k   and 


(6.43) 


k   are  known  as  the 


The  maximum  transmittance  to  linearly  polarized 
:   and  the  minimum  transmit tanc 
90°)  is   k  .   The  relationships  (6.43)  imply  that 


light  is   k   and  the  minimum  transmittance  (when  the  direction  of  polarization  is  rotated 


q2  =  s2  -  d2  -  p2. 


(6.44) 


We  use  the  word  "filter"  here  and  in  the  rest  of  this  chapter  in  a  general  sense  to  refer 
to  any  optical  element  which  may  modify  the  distribution  of  radiation  with  respect  to  any 
parameter,  including  polarization. 

2Equation  (6.42)  can  also  be  used  to  describe  the  polarizing  effects  of  reflecting  surfaces 

[4,14].   If  the  coordinate  system  described  in  the  footnote  on  page  9  is  used  then 

,  o 


3Nominal  values  of   k..   and   k„   have  been  published  [2]  for  the  common  Polaroid  polarizers, 

30 


Thus  we  really  have  only  three  parameters  to  measure  —  say   s,  d,  and  p  —  in  order  to 
define  the  Mueller  matrix  p  completely.   Note  that  the  matrix  includes,  as  special  cases, 
both  perfect  polarizers   (k  =  1,   k  =  0)   and  perfect  retarders   (k  =1,   k  =  l)  .   It 
also  provides  for  the  most  important  imperfections  encountered  in  real  linear  polarizers 
and  retarders. 

We  introduce  the  notation  F(<j>)   to  indicate  the  result  of  a  measurement  with  our 

polarization- indifferent  detector  of  the  transmittance  to  unpolarized  radiation  of  filter 

F  when  it  is  turned  to  an  orientation  angle   d>   on  its  mount.   Likewise   GF(d>  ,*_)   is  the 

g   f 
measured  transmittance  of  the  pair  of  filters   F   and   G  at  angular  settings  of   <j>   and 

<()  ,  respectively.   Our  notation  also  indicates  the  direction  of  radiation  travel  —  from 

right  to  left  through  the  symbols;  i.e.  through  filter   F   first,  then  through  filter   G 

to  the  detector. 

We  first  measure  the  transmittances   F(d>_),  G  (<i>  ),   and  H(<i,  )   of  the  three  individual 

r      g  h 

filters.   From  eq.  (6.42)  we  see  that  these  are  given  by 


F(+f) 


sf,   G 


)  =  s  ,   and  H  (<j>,  )  =  s, 
;     s  h     h 


where  the  subscripts  on   s   identify  the  filter  to  which  it  pertains.   This  result  should 

be  independent  of  i>    ,    i>    ,    and  d>,  ,   so,  in  principle,  a  measurement  at  any  orientation 
f   g       h 

should  suffice.   However,  a  real  filter  will  probably  have  small  inhomogeneities  which  may 
move  into  and  out  of  the  beam  as  the  filter  rotates.   In  this  case  the  best  that  can  be 
done  is  to  obtain  an  average  value  of   s  by  sampling  many  orientations   <£.   These  averages 
of   F((|>  )  ,  etc.,  we  will  simply  write  as   F,  G,  and  H.   Thus, 


s_  =  F,   s   =  G,   and   s,  =  H 

f        g  h 


(6.45) 


We  now  measure   the   transmittance     FH(<J>f,<j>    )     of   a  pair   of    the   filters.      From  eq.    (6.42) 


we  have 


F(*f)*H<*h>  = 


[sf*sh  +  Wcfh] 
LVVcf  +  VVVcfh  +  VVVsfh] 

[Wsf  +  sf-VVcfh  ~  Pf-VVsfh] 
[Wsfh] 


(6.46) 


where  we  have  shown  only  the  first  column  of  the  product  transmittance  matrix  because  it 
turns  out  that  this  is  all  that  enters  into  the  final  expressions.   We  have  also  used  the 
shorthand  notation 


Cfh  =  cos2[(cDf  -  d,f°)  -  Uh  -  c)>h0)] 


and 


S,,  =  sin2[Uf  -  4>°)    -  ( 


fh 


f    Yf 


Kh   Th 


°>] 


31 


The  transmittance   FH(d),.,<f),  )  is  just  the  (0,0)   element  of  this  matrix  or: 

t   n 


ra<*f.*h>  =  Sf'Sh  +  WCfh 


By  expanding   C,.,   we  can  rewrite  this  as 
1  h 


FH< 


>fl<j>h)    =   FHX  +  FH2-cos2(<t>f 


|>h)    +  FH3.sin2(4>f   -   <j>h) 


(6.47) 


where 


FH1  "  Sf'Sh 


FH2  =  df-dh'cos2(*f°  -  4»h°) 


and 


FH3  =  df.dh.sin2(*f°  -*h°) 


If  we  measure   FH  ($  ,<j>  )   at  two  or  more  sets  of  angles   <J>   and   <f>   we  will  be  able  to 
solve  the  resulting  set  of  equations  (6.47)  —  either  by  the  usual  methods  for  simultaneous 
equations  or,  better,  by  least  squares  fitting  —  for  the  coefficients   FH   and   FH  . 
FH    of  course,   is  already  known  from  eq.  (6.45).1   In  this  way  the  orientation- independ- 
ent coefficients   FH   and   FH   are  measured.   From  these  we  obtain 


tan2(< 


^h  >  =  FH3/FH2 


and 


V 


Vdh=:  VFH22  +FH32 


(6.48) 
(6.49) 


By  repeating  this  process  for  the  other  two  possible  ways  of  pairing  the  three  filters  we 
obtain  the  other  polarization-axis  angle  differences  and  the  other   d*d   products  like 
eq.  (6.49).   The  three  equations  of  the  form  (6.49)  can  be  further  solved  for  the  individ- 
ual values  of  d: 

-i  u 


df  = 


d   = 


dh  = 


2     2\  / 2     —  2 

GH2   +  GH3j-(^FH2   +  FH3 

—  2    —  2 
HG2  +.HG3 


2    —  2 
HG2   +  HG3 


GF22  +  GF32 


•  2     —  2 
GF2   +GF3 


FH22  +  FH32 
FH22     +  FH32J-(hG22     +  HG32 


(6.50) 


If 


FH   is  determined  from  eqs.  (6.47)  and  does  not  agree  with  F»H   [eq.  (6.45)]  then 


the  model  of  eq.  (6.42)  may  be  inadequate.   Most  likely  some  of  the   (0,3)   or 
elements  don't  vanish  as  assumed  in  eq.  (6.42). 
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(3,0) 


So  far  our  measurements  have  given  us  the  orientation  of  the  three  filter  polarization 
axis  directions   <J>   with  respect  to  one  another  but  have  not  revealed  the  absolute  orienta- 
tion of  any  one  of  them.   This  is  obviously  because  there  is  no  unique  transverse  reference 
direction  associated  with  our  unpolarized  measurement  beam.   The  absolute  polarization  axis 
orientation  can  be  determined,  if  desired,  by  passing  a  beam  known  to  be  partially  hori- 
zontally linearly  polarized  through  one  of  the  filters  and  rotating  the  filter  for  maximum 
transmission.   At  this  orientation   <}>  =  <J>   which  then  serves  to  determine  the  <\>        values 
of  the  other  two  filters  from  eqs.  (6.48).   A  suitable  source  of  partially  horizontally 
polarized  light  is  light  reflected  from  a  glass  surface  oriented  so  that  the  plane  of 
incidence  (the  plane  defined  by  the  incident  and  reflected  rays)  is  vertical. 

The  final  measurements  are  of  the  three-filter  transmittances   HGF(<j>  ,<J>  ,$>f).      From 
eqs.  (6.46)  and  (6.42)  we  find,  after  some  simplification,  that  the   (0,0)   element  of  the 
three-filter  product  matrix  is: 


HGF(<j>,  ,(f>    ,((>,,)    =    s    -s    >sc   +  s    -d    -d    -C    c   +  d-d    -s    -Q 

h     g     f  hgf  hgfgf  hgfhg 


+  d    «s    -d    «C   C'C,       -  d    «p    «d    «S    C'S,     . 
h     g      f      gf      hg  h   rg      f      gf      hg 


(6.51) 


As  before  C   =  cos2  [(<J>   -  cf>  °)  -  (<j>.  -  <j>  °)]  =  cos2  [(<f>   -  $  )  -  (<J>  °  -  <j>  °)]  ,  etc.   Every- 
thing in  this  equation  is  now  known  [from  eqs.  (6.45),  (6.48),  and  (6.49)]  except   p  ,   so 
p   can  be  calculated  directly: 


P„  = 


+  G 


H-G-F  -  HGF(<h  ,<(>  ,«j)  )  +  H-(GF„2  +  GF2)2-C   +  F-(HG2  +  HG.2)^-C, 
h   g   i I 3 g_f 2 J     hg 

2    —  2  k 

(FH„  +  FH_  )   -S   -S 
2     3     gf  hg 

gf   hg 

S  f'SU        ' 

gf   hg 


(6.52) 


Again  this  result  should  be  averaged  over  many  settings  of   4>  ,  (f>  ,  and   (j>    in  order  to 
minimize  errors  due  to  filter  inhomogeneities.   Similar  three-filter  measurements  with 
filters   F  and  H   in  the  middle  will  permit   pf   and   p   to  be  determined.   Finally  the 
values  of   q   can  be  calculated1  from  eq.  (6.44)  and  with  that  the  three  matrices  corre- 
sponding to  eq.  (6.42)  are  completely  evaluated. 


The  algebraic  sign  of   q   is  not  determined  by  the  procedure  we  have  outlined  here.   In 
many  applications  this  will  be  of  no  consequence.   The  ambiguity  is,  in  part,  due  to  the 
arbitrariness  in  the  definition  of  circular  handedness.   The  sign  of   q   can  be  determined 
if  a  beam  containing  circularly  polarized  light  of  specified  handedness  is  passed  through 
the  filter  and  the  output  is  analyzed  with  a  linear  polarizer.   See  also  the  discussion  of 
these  sign  ambiguities  in  Appendix  7. 
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The  procedure  we  have  outlined  applies  only  to  filters  which  are  at  least  partial 
linear  polarizers.   This  is  because  the  right-hand  side  of  eq.  (6.48)  can  become  indeter- 
minate  (0/0)   for  a  pure  retarder  or  a  simple  attenuating  spectral  filter.   Any  such 
component  can,  however,  be  measured  by  putting  it  in  the  middle  of  a  three-filter  train 
between  two  polarizing  filters  which  have  already  been  characterized.   Equation  (6.51)  is 
still  valid  and  the  symbol   G  now  refers  to  the  unknown  filter.   The  quantity   s   can  be 

obtained  by  the  single  filter  transmittance  measurement  as  before.   This  leaves   d  ,   p  , 

o  o 

and   <j>     (appearing  in  C  f,   etc.)  to  be  determined.   The  angle   $    can  be  established 

g  o  g   o  8 

by  setting   <j>   =  <j>,  +  <j>f  -  <j>    (i-e.,  the  polarization  axes  of   F   and  H  are  made 

parallel  to  each  other)  and  rotating  the  unknown  filter  G  until  a  maximum  transmittance 
is  obtained.  With  the  polarization  axes  of  F  and  H  parallel  to  each  other  the  trans- 
mittance is 

HGF(<f>.  ,4>  ,<fr.  +  <J>  °  -  <f>,°)  =  s,  -s  -s^  +  d,  -p  -d,  +  (s  -d.  +  d,  •■s_)  »d  -C. 

Yh'Yg  Th    Tf    Th     h  g   f    h  rg   f     h   f    h   f    g  hg 


+  d,  'd.-(s   -  p  )-C  2 
h   f   g    g   hg 


(6.53) 


This  expression  has  a  maximum  at   d>   =  <f>,  +   d>    -  6,  .  1   This  serves  to  determine  the  angle 

g    h    g     h 
(j)  ,  which  once  again  should  be  averaged  for  many  settings  of   d>,  .   Once   d>    is  known  the 
g  h  g 

other  parameters  of  filter  G  can  easily  be  obtained  by  making  many  measurements  and  using 

linear-least-squares  fitting  in  eq.  (6.51)  or  by  making  measurements  at  two  suitable  angles, 

4>  ,   and  solving  the  resulting  two  equations  given  by  eq .  (6.53)  for  the  unknowns   d   and 

p  .   q   then  follows  as  before  from  eq.  (6.44).   If  the  transmittance,  eq.  (6.53), 

exhibits  no  significant  dependence  on  <t>,     -   <t>        after  averaging  over   d>,   then  d   =0, 

h    g  h         g 

and   p   =  s  .   Such  a  filter  is  just  an  attenuating  spectral  filter. 

Needless  to  say,  with  real  measurements  on  real  filters  some  uncertainties  in 
determining   s,   d,   p,   and   q   can  be  expected.   For  example,  if  1%  measurements  are  being 
made  of  a  good  linear  polarizer  the  uncertainties  in   s,   d,   p,   and   q  will  all  be 
similar  and  of  the  order  of  s/100.   But  since   p  ^  q  ^  0   the  relative  uncertainty  in  these 
two  parameters  may  appear  alarmingly  large.   Fortunately,  the  relative  contribution  to  any 
final  instrument  response  from  a  single  matrix  element   T . .   of  an  optical  component  cannot 

be  larger2  than  the  fraction   |t../T„_|.   Thus,  in  the  example  of  the  1%  measurements  of  a 

°  '  ij   00 ' 

linear  polarizer  the  uncertainties  in  p   and   q   can  never  contribute  more  than  about  a 
1%  uncertainty  to  any  measurement  made  using  this  polarizer. 


If   d   is  small  but  not  zero  there  may  be  a  lesser  maximum  90°  removed  from  this  angle. 


2This  is  true  if  the  ultimate  instrument  detector  responds  to  a  property  of  the  electro- 
magnetic radiation  such  as  power,  linear  momentum,  or  photon  number,  but  need  not  be  true 
if  it  responds  to  a  property  characteristic  of  a  particular  polarization  state  such  as 
angular  momentum. 
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Once  the  Mueller  matrix  elements  have  been  determined  for  two  real  linear  polarizers 
and  two  approximate  quarter-wave  plates  the  polarizer  and  wave-plate  pairs  can  be  used  to 
generate  and  analyze  a  wide  range  of  polarization  states  of  light.   And  this  in  turn  makes 
possible  the  characterization  of  more  complicated  Mueller  matrices  —  in  particular  those 
of  circular  polarizers,  which  are  not  included  in  our  general  imperfect  filter  expression 
p   in  eq.  (6.42)  because  we  have  set   F   =  F  .  =0.   We  thus  see  how  we  can  start  with 
real,  inexpensive,  polarizing  filters  and  build  up  a  complete  polarization  measurement 
capability. 

POLARI ZER- ATTENUATORS 

A  three-polarizer  train  is  sometimes  used  as  an  absolute  or  predictable  variable  opti- 
cal attenuator  [15,16].   The  first  and  last  polarizers  of  this  train  are  aligned  with  their 
polarization  axes  parallel  to  each  other  and  the  middle  element  is  rotated  to  achieve  the 
desired  degree  of  attenuation.   Equation  (6.53)  describes  just  this  arrangement.   The  trans- 
mittance  of  the  train  relative  to  its  maximum  transmittance  is  therefore: 


HGF(cJ>°,<j>      <t>°) 

t    =  — - — S  =    [S    .s    .-     +  d,  -d.-p     +   (s,-d£  +  d    -s,).d    -C, 

o,o,  ox  Lhgf  h      f    *g  hf  hfghj 


HGF((|)h,(J.g,<j)f  ) 


+  d,  -d.-(s      -  p   )'C,    2]  /[s,  -s    -s^  +  s    -d    -d^  +  d    -s    -d^ 
hfvg        Vhghgf  hgf  hgf 


(6.54) 


+  d-d  -sJ 
h  g  f 


his  expression  for  the  relative  transmittance,  of  course,  is  valid  only  for  unpolarized 
light  and  a  polarization- indifferent  detector.   If  the  polarizers  are  all  perfect  linear 

P, 


polarizers  we  have   s^  =  s   =  s,  =  1,   d,.  =  d   =  d,  =  1,   and   p,. 

f    g    h        f    g    h  vl 


=  d,  =0.   Then 
'  h 


1  +  2C  +  C 


=  cos  (c 


')  • 


If  the  outer  polarizers  are  perfect  and  the  middle  component  is  a  perfect  half-wave  plate 
instead  of  a  polarizer  the  only  change  is  that   d   =0   and  p   =  -1.   This  results  in 


2C 


-£-  =  cos22(<J)   -  <j>  °) 


In  both  these  cases  the  light  emerging  from  the  attenuator  is  always  completely  linearly 

polarized  in  the  direction  determined  by  the  final  polarizer.   The  relative  transmittance 

depends  only  on  the  easily  measured  angle   <J>  . 

g 

Unfortunately,  with  real  optical  elements  the  values  of   s,   d,   and  p  must  be  known 
if  the  transmittance  is  to  be  calculated  accurately  as  a  function  of   <j>    from  eq.  (6.54). 
Moreover,  the  state  of  polarization  of  radiation  emerging  from  the  attenuator  will  depend 
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weakly  upon  the  angle  $        and  upon  the  state  of  polarization  of  the  entering  radiation. 
Thus,  the  Mueller  matrix  of  the  rest  of  the  optical  system  must  be  known  and  additional 
matrix  elements  of  the  attenuator  must  be  evaluated  if  the  effective  relative  transmittance 
is  to  be  calculated  exactly. 

DEPOLARIZERS 

We  briefly  mentioned  depolarized  radiation  in  the  early  part  of  this  chapter  and 
would  now  like  to  say  a  few  words  about  the  means  of  achieving  this.   The  Mueller  trans- 
mittance matrix  for  a  depolarizer  is  obviously  proportional  to 


10  0  0 

0   0  0  0 

0   0  0  0  (6>5M 

0   0  0  0 


Unfortunately  no  single  optical  element  exhibits  this  property.   However,  if  we  consider  a 
complete  measurement  situation  with  a  source  radiance  vector   |_,   an  optical  component  with 
Mueller  matrix  [),   and  a  radiometer  with  responsivity  matrix  ^P  ,  with  the  elements 
(responsivity  coefficients)   R^j   Rni ,   etc.,  we  have  for  the  radiometer  output  signal: 


S  = 


At 


AX 


AA 


Aw 


/^•D-L,aw'cose'dA'dX'dt 


(6.56) 


where   Aio,   AA,   AA,   and   At   are  the  resolution  intervals  of  the  radiometer  in  these 
parameters.   (As  in  our  discussion  of  the  calibration  of  a  radiometer,  we  actually  use  only 
the  upper  row  of  the  responsivity  matrix  /^    to  obtain  the  output  signal   S.)   If  the 
responsivity  of  the  radiometer  and  the  Stokes  spectral  radiance  vectors  of  the  source  are 
sufficiently  uniform  over  the  radiometer  resolution  intervals,  they  can  be  removed  from 
this  integral  to  give 


S  - 


£• 


At 


AA 


AA 


D'du'cosfi 'dA'dA-dt 


Aco 


AfAA'AA'AO 

=  /R-  <D>#L'At'AA'AA'A0 


•L'Af  AA-AA'AG 


(6.57) 


where 


<D>  = 


<D00>  ^Ol*  <D02>  <D03> 

<D1Q>  <DU>  <D12>  <D13> 

<D2Q>  <Dn>  <D22>  <D23> 

<D3()>  <D31>  <D32>  <D33> 
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and   <D  .  .  >  = — — — — — 

lj     At-AA'AA'AG 


At 


AAJ  AA 


.'■.w 


D.  •dorcos6'dA«dA'dt.   Now  we  see  that  if  all   <D . . > 


except   <Dnn>   can  be  made  to  vanish,  the  Mueller  matrix  <Q>  will  become 


00 


<D> 


<D      > 
00 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

(6.58) 


which  is  the  form  for  a  depolarizer.   Most  depolarizers  consist  of  two  retarders  [eq.  (6.20)] 
placed  one  behind  the  other.   The   (0,0)   element  of  the  resulting  product  matrix  is   1   and 
all  the  other  elements  are  zero  or  are  trigonometric  functions  of  the   <£ ' s   and   6's   of  the 
two  retarders.   By  carefully  choosing  the   <j>   values  (or  the   6   values)  and  averaging  over 
a  range  of   6's   (or   <j>'s)   these  trigonometric  matrix  elements  can  be  made  negligibly 
small.   The   6's   (or   <j>'s),  over  which  the  device  averages,  are  arranged  to  vary  in  time, 
wavelength,  position  or  direction.   Depolarizers  which  depolarize  by  averaging  over  wave- 
length (^100  nm)  [17]  or  over  time  (^1  sec)  [17]  or  over  the  cross-sectional  area  of  the 
beam  (^1   cm2)  [3]  have  been  used  and  will  be  described  more  completely  in  Part  II.   Depolar- 
izers are  frequently  known  as  pseudo-depolarizers  or  as  polarization  scramblers. 

A  SIMPLIFIED  CHARACTERIZATION  of  POLARIZATION 

Many  radiometric  measurements  involving  polarized  light  are  carried  out  using  just  two 
measurements  —  one  with  a  linear  polarizer  set  in  a  horizontal  orientation  in  front  of  the 
detector  and  the  second  with  the  polarizer  turned  90°.   Let  us  see  what  the  limitations 
of  such  measurements  are  and  under  what  conditions  they  will  be  sufficient  to  characterize 
a  light  ray.   In  order  not  to  complicate  matters  unduly,  let  us  assume  that  the  polarizer 
is  perfect  so  that  its  Mueller  transmittance  matrix  in  the  horizontal  orientation  is: 

110  0 

110  0 

0   0   0  0 

0   0   0  0 


P(0°)  = 


and  in  the  vertical  orientation: 


P(90°)  = 


1-10  0 

-110  0 

0   0   0  0 

0   0   0  0 


We  will  also  assume  that  the  detector  is  polarization  indifferent  or  has  been  compensated 
or  calibrated  for  horizontally  and  vertically  polarized  radiation.   When  an  arbitrarily 
polarized  ray  with  Stokes  vector 
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passes  through  the  linear  polarizer  in  the  two  orientations,  we  have 


L*(o°)  =  P(o°).[_  = 


Lo  +  Li 
Lo  +  Li 

0 
0 


and 


L'(90°)  =  P(90°)-L  = 


Lo"Li 
"Lo  +  Li 

0 
0 


The  detector  will  record  signals  proportional  to  L  '(0°)  =  L_  +  L    and   L  ' (90°)  =  L~  -  L,  , 
respectively.   Thus  the  original  components   L   and   L   can  be  calculated  from 


%[L'(0°)  +  L  '  (90°)]   and  %[L  '  (0°)  -  L  '  (90°)].   If  we  can  assume  that   L 


0   then 


the  light  beam  has  been  completely  characterized.   Usually  L  ,  the  spectral  radiance,  and 
the  degree  of  polarization  [eq.  (6.6)] 


P  = 


L0'(0°)  -  L0'(90°) 
LQ     L0'(0°)  +  L0'(90°) 


are  the  quantities  reported.   Thus  we  see  that  the  two-measurement  description  of  polarized 
light  is  limited  to  those  situations  in  which  (as  we  could  have  guessed)  there  is  no  45° 
(or  135°)  linearly  polarized  component  and  no  circularly  polarized  component.   The  elimina- 
tion of  the  +45°-preference  component  can  always  be  arranged  by  suitable  choice  of  the 
reference  "horizontal"  direction.   Note,  however,  that  unless  one  has  some  a  priori   knowledge 
of  this  direction  this  entails  the  equivalent  of  at  least  one  additional  measurement  to  find 
this  direction  and  there  is  thus  no  real  saving  in  attempting  to  eliminate  this  component. 
In  many  situations,  such  as  scattering  studies  of  isotropic  samples  irradiated,  first,  by 
radiation  polarized  in  the  plane  of  incidence  and  then  by  radiation  polarized  perpendicular 
to  the  plane  of  incidence,  a  reference  direction  parallel  to  or  perpendicular  to  the  plane 
of  incidence  can  be  taken  and  the  +45°-preference  component  may  then  be  safely  neglected 
(but  only  if  the  sample  is  isotropic) .   The  neglect  of  the  circularly  polarized  component 
is  probably  usually  justified  in  most  radiometric  measurements  where  uncertainties  of  a 
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percent  or  so  are  tolerable.   For  the  highest  accuracies  circular  polarization  must  be 
included,  however,  because  it  can  creep  in  whenever  partially  polarized  light  at  non-normal 
incidence  is  reflected  from  mirrors  or  transmitted  through  films.   Thus,  any  but  the 
simplest  instruments  will  be  likely  to  exhibit  some  circular-polarization  selectivity. 

When  circular  polarization  can  be  neglected  most  of  the  discussions  of  this  chapter 
can  be  greatly  simplified.   For  example,  no  measurements  involving  retarders  such  as  quar- 
ter-wave plates  need  be  considered;  sources,  optical  components  and  radiometers  can  be 
completely  characterized  using  only  linear  polarizers.   Also  the  procedure  outlined  in 
eqs.  (6.42)  to  (6.53)  for  measuring  the  Mueller  matrix  of  a  polarizer  can  be  simplified  by 
omitting  the  three-filter  measurements  because  we  can  now  set   q  =  0   and   p2  =  s2  -  d2 
for  each  filter. 

EXPERIMENTAL  CONCERNS 

We  have  tried  in  this  chapter  to  present  a  complete,  exact  treatment  of  polarization 
and  an  indication  of  how  this  treatment  can  be  applied  in  the  laboratory.   We  have  not 
mentioned  the  experimental  problems  which  make  high  accuracy  polarization  measurements 
difficult  —  these  properly  belong  in  later  chapters  where  we  deal  with  specific  radiometric 
instruments.   However,  some  brief  mention  of  some  of  these  problems  seems  appropriate  here. 
The  proper  experimental  inclusion  of  polarization  in  measuring  instrumentation,  as  we  have 
seen,  requires  the  insertion  into  the  optical  path  of  one,  two,  or  three  polarizing  filters 
of  one  kind  or  another.   These  must  be  rotatable  and  it  must  be  possible  to  measure  the 
angles  of  rotation,  or  at  least  to  set  each  filter  reproducibly  at  several  preselected 
angular  positions.   Unless  the  filters  are  absolutely  uniform  there  will  be  a  dependence  of 
filter  transmittance  on  angular  orientation  which  may  be  unrelated  to  the  polarization 
property  being  investigated.   For  anything  more  than  qualitative  observations  it  is  there- 
fore important  that  polarization  measurements  be  carried  out  and  suitably  averaged  over 
many  orientations  of  the  polarizing  components.   Another  experimental  difficulty  common  to 
all  measurements  using  multiple  filters  in  tandem  but  which  is  exacerbated  with  polarizing 
filters  is  that  of  filter-filter  interref lections .   Such  interref lections  have  the  effect 
of  increasing  the  apparent  transmittance  of  the  filter  train  and  thereby  introducing  errors. 
In  the  case  of  polarizing  filters  this  problem  can  be  eliminated  by  the  usual  technique  of 
tilting  the  filters,  but  the  axis  of  rotation  of  the  filters  must  remain  parallel  to  the 
beam  axis.   Thus,  the  filter  must  wobble  as  it  is  rotated.   If  a  tilted  polarizing  filter 
were  rotated  about  an  axis  perpendicular  to  its  face,  then  entering  radiation  would  be 
partially  linearly  polarized  at  a  fixed  orientation  by  the  tilted  first  surface,  resulting 
in  observations  of  partial  polarization  even  of  unpolarized  light.   If  the  tilt  rotates 
with  the  filter,  however,  then  this  polarizing  effect  becomes  simply  an  integral  property 
of  the  filter  and  will  be  accounted  for  in  the  measured  Mueller  transmittance  matrix  ele- 
ments.  Of  course,  the  filter  must  always  be  used  at  the  same  tilt  orientation.   Another 
complication  which  may  require  experimental  attention  is  that  the  lateral  displacement  of 
the  beam  as  it  passes  through  the  tilted  filter  now  rotates  with  the  filter. 
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SUMMARY  of  CHAPTER  6. 

Polarization  can  be  included  rigorously  in  classical  radiometry  (geometrical  optics) 
by  treating  spectral  radiance  as  a  four-component  mathematical  column  vector  called  a 
Stokes  vector 


Lx- 


X,0 

Vi 

LA,2 

Lla,3 


(6.8a) 


The  four  components  of   |_   are  known  as  Stokes  components.   The  Stokes  vector  and  its 

components,  in  general,  depend  upon  ray  position  and  direction  and  upon  wavelength,  i.e., 

upon  x,y,0,(j>,  and  A.   The  first  Stokes  component   L   n   is  the  spectral  radiance  which  a 

A  ,  (J 

polarization- indifferent  detector  would  measure.   It  is  the  measure  of  the  spectral  radiant 
flux  associated  with  the  ray  and  is  identical  to  the  spectral  radiance  discussed  in  earlier 

chapters.   Each  of  the  other  three  components  of   |_   represents  the  excess  of  one  form  of 

A 

polarization  over  its  complementary  form.   L   .,   describes  the  excess  of  horizontal  over 

A,i 

vertical  linear  polarization,   L   „   describes  the  excess  of  +45°  over  +135°  linear  polar- 

A,Z 

ization  and  L   _   describes  the  excess  of  right  over  left  circular  polarization.   The 

A,  J 

spectral  radiance  L   n,   is  the  sum  of  two  parts,  a  polarized  part,   L    ,   given  by 

A  ,  U  A  ,p 


,p    1     A,. 


+  L 


A, 2 


+  L 


A, 3 


(6.5) 


and  an  unpolarized  part 


The  fraction 


A ,u    A ,0    A ,p ' 


A,p   A,0 


(6.6) 


is  known  as  the  degree  of  polarization. 

In  this  description  of  polarized  radiation  the  transmittance  or  propagance  of  an 
optical  path  is  represented  by  a  4  x  4  matrix  known  as  a  Mueller  matrix.  For  example 
(compare  eq.  (3.14)  [1]) 


l»  -  t-l, 


(6.9a) 


relates  the  spectral  radiance   |_   at  one  point  on  a  ray  and  the  propagance   )   along  the 

A 

path  between  that  point  and  a  second  point  to  the  spectral  radiance   |_  '   at  the  second 

A 

point.   The  Mueller  propagance  matrix     in  this  equation  is  given  by 
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so  that  eq.  (6.9a)  is  equivalent  to 


LJ 


T00'L0  +  T01-Ll  +  T02'L2  +  T03'L3 


Ll  "  T10-L0  +  T11'L1  +  T12-L2  +  T13'L3 


L2   =  T20'L0  +  T2l'Ll  +  T22*L2  +  T23*L3 


(6.7) 


L3    T30'L0  +  T3l'Ll  +  T32'L2  +  T33*L3 

where  we  have  dropped  the  subscript   A   for  simplicity. 

When  polarization  is  included,  the  spectral  radiance  of  a  ray  of  radiation  transmitted 
by  a  train  of  optical  components  depends  upon  the  order  in  which  the  propagating  flux  in 
the  ray  encounters  the  components.   This  is  reflected  in  matrix  equations  by  the  order  of 
writing  the  Mueller  matrices  and  the  Stokes  vectors.   The  entering  Stokes  vector  is  written 
on  the  right  and  the  Mueller  matrices  are  written  from  right  to  left  in  the  order  in  which 
the  optical  components  they  represent  are  encountered.   The  resulting  matrix  equations  can 
then  be  evaluated  by  the  usual  matrix  multiplication  rules  (Appendix  6) . 

The  Mueller  transmittance  matrix  of  an  ideal  linear  polarizer  is  given  by 


P(4>)  =  h 


1 

c 

s 

0 

c 

c2 

s-c 
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0 
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0 

(6.15a) 


where   C  =  cos2<J>,   S  =  sin2<{>,   and   <j>   is  the  angle  between  the  polarization  axis  of  the 
polarizer  and  the  horizontal  polarization  reference  axis  of  the  ray.   When  a  ray  of  arbi- 
trary spectral  radiance  is  passed  through  such  a  polarizer  the  polarization  state  of  the 
transmitted  ray  can  be  predicted  by  performing  the  matrix  multiplication  indicated  in 
eq.  (6.9a): 


|_i  =  P(<|>).|_  =  P(< 


_  3_ 


k(.LQ    +    L1-C  +  L2«S) 


(6.17) 
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The  transmitted  ray  has  a  spectral  radiance  of  ^(L.  +  L  -cos2<}>  +  L  «sin2<J>).   It  has  no  net 
circularly  polarized  component   (L  '  =  0)   but  is  completely  linearly  polarized  at  the  angle 
<j>,  [P  =  1,  from  eq.  (6.6)]  . 

The  Mueller  transmittance  matrix  of  an  ideal  linear  retarder  is 


Res.*)  = 


1 

0 

0 

0 

0 

C2  +  S2-cos6 

S-C- (1  -  cos6) 

-S-sin6 

0 

S-C-  (1  -  cos6) 

S2  +  C2-cos6 

C-sin6 

0 

S-sin6 

-C-sin6 

cos6 

(6.20) 


where   C  =  cos2cj>,   S  =  sin2cf>,   and   <f>   is  the  angle  between  the  fast  axis  of  the  retarder 
and  the  horizontal  reference  axis  of  the  ray.   6   is  the  retardance  of  the  device  and  is 
proportional  to  the  thickness  of  the  retarder  plate.   If   6  =  it/2  [rad]   the  retarder  is 
called  a  quarter-wave  plate,  which  has  a  Mueller  transmittance  matrix 


QUO  = 
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c2 

S-C 
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s2 

c 

0 
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(6.21) 


The  combination  of  an  ideal  quarter-wave  plate  and  an  ideal  linear  polarizer  can  be 
used  either  to  determine  the  state  of  polarization  of  an  arbitrarily  polarized  ray  or  to 
produce  any  given  state  of  polarization.   If  the  radiation  passes  first  through  a  rotatable 
quarter  wave  plate  at  angle   a   and  then  through  a  fixed  polarizer  at  angle   3   the  trans- 
mitted ray  is  given  by 


L'  =  P(3)-Q(a; 


Jo 


LL3J 


=  h[L     +   L  • cos2a«cos2(8-a)  +  L  •sin2a«cos2($-a)  +  L  •sin2($-a)] 


(6.23) 


1 
cos23 
sin2g 

0 


If  the  spectral  radiance  of  the  transmitted  ray  is  measured  at  four  suitable  orientations 
a   of  the  quarter-wave  plate  the  four  simultaneous  linear  equations  of  the  form  of 
eq.  (6.23)  can  be  solved  for  the  Stokes  components  L  ,  L..  ,  L„ ,  and  L„   of  the  original 
ray.   Since  the  transmitted  beam  is  always  completely  linearly  polarized  at  the  fixed 
orientation   3,   the  measurement  can,  in  principle,  be  carried  out  by  a  detector  of  any 
arbitrary  polarization  selectivity.   For  absolute  measurements,  however,  a  calibration  of 
the  detector  will  be  required  using  (usually)  a  differently  polarized  source.   In  this  case 
it  will  be  simpler  if  the  detector  is  polarization-indifferent  so  that  the  spectral-radiance, 
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Ln' ,   is  the  quantity  measured, 


If  the  radiation  passes  first  through  the  polarizer  and  then  through  the  quarter-wave 
plate,  the  transmitted  ray  will  be  described  by 


|_'    =   Q(a)  -P(3) 


h(LQ  +  L1-cos23  +  L2-sin2g) 


cos2a*cos2(a-$) 

sin2cccos2(a-3) 

sin2(a-3) 


(6.35) 


By  suitable  choice  of  the  orientation  angles  of  the  quarter -wave  plate  and  polarizer  any 
pure  polarization  state  from  circular   (|a-g|  =  45°)   to  linear   (a-3  =  0°)   can  be  realized, 

The  final  electrical,  thermal,  or  chemical  response,   S,   of  a  radiometer  depends  upon 
all  four  Stokes  components  of  the  entering  radiance: 


S  = 


AAJ  AA 


Aw 


(RQ0-L0  +  Rq-^I^  +  Rq2-L2  +  R03-L3)-dco-cos6-dA-dA, 


(6.36) 


The  quantities   AA,  AA,  and  Aoj  are  the  acceptance  intervals  of  the  radiometer.   If  these 
acceptance  intervals  are  sufficiently  small  that  the  Stokes  spectral-radiance  components 
can  be  assumed  constant  and  removed  from  the  integrals  in  eq.  (6.36)  then  we  can  write 


S=*00'L0+/?01-Ll+i?02-L2+/?03-L3- 


(6.36a) 


The  values  of  the  responsivity  factors  R  can  be  determined  experimentally  from  observa- 

tions of  the  radiometer  output  for  a  beam  in  four  different,  known  states  of  polarization 
such  as  can  be  obtained  from  the  polarizer  and  quarter-wave  plate  combination  described  by 
eq.  (6.35).   If  the  first  optical  component  of  the  radiometer  is  a  perfect  diffuser, 
eq.  (6.36)  reduces  to 


S  =  R 


where 


00  0 


AA 


AAJ  Ago 


(6.36c) 


L  'dwcose  *dA*dA   is  the  total  radiant  flux  falling  on  the  diffuser 


within  the  wavelength  band  pass  of  the  radiometer.   If  neither  of  these  simplifications  is 
appropriate  the  dependence  of  the  coefficients   R    upon  ray  position,  direction  and  wave- 
length can,  in  principle,  be  determined  by  probing  the  radiometer  with  rays  of  known  polar- 
ization state.   With  these  functions  known,  then,  a  radiometer  measurement  of  an  unknown 
beam  provides  one  relationship  among  integrals  of  the  four   L.(x,y,0,<j>,A)   functions.   If 
the  use  of  the  radiometer  is  restricted  to  comparisons  of  beams  with  identical  relative 
spectral-radiance  distributions  then  such  detailed  knowledge  of  its  responsivity  is  unnec- 
essary and  four  measurements  of  different  known  polarization  states  will  characterize  it 
for  this  purpose. 
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The  ideal  Mueller  transmlttance  matrices  of  eqs .  (6.15)  and  (6.20)  are  only  approxi- 
mations to  the  transmlttance  matrices  of  real  optical  components.   Many  real  linear 
polarizers  and  retarders  can,  however,  be  adequately  represented  by  Mueller  matrices  of  the 
form 


F(+)  = 


s-C2  +  p-S2 


s        d-C 
d-C 

d-S  (s-p)-S-C 

0        q-S 


d-S 

(s-p)-S-C 

s-S2  +  p-C2 

-q-C 


0 
-q-S 
q-C 

P 


(6.42) 


where 


and 


C  =  cos2(<})-cf>0)  ,   S  =  sin2(<j>-c|>0) 


q2  =  s2  -  d2  -  p2. 


(6.44) 


The  values  of   s,  d,  and  p   as  well  as  the  filter  polarization  axis  direction  $        can  be 
determined  experimentally  by  appropriate  transmittance  measurements.   If  no  sources  of 
known  variable  polarization  state  are  available  it  is  possible  to  measure  these  parameters 
with  sufficient  accuracy  for  most  radiometric  applications  by  making  transmittance  measure- 
ments in  unpolarized  light  on  three  or  more  such  objects  as  a  function  of  their  orientation 
angles.   Once  the  Mueller  transmittance  matrices  of  a  polarizer  and  an  approximate  quarter- 
wave  plate  have  been  established  they  can  be  used  for  analyzing  or  generating  a  wide  range 
of  arbitrary  polarization  states  by  generalizing  eqs.  (6.23)  and  (6.35). 
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Appendix  6.    Matrix  Multiplication 

by  John  B.  Shumaker 
In  this  appendix  we  review  the  rules  for  matrix  multiplication.   We  will  confine  our 
attention  to  the  4-component  Stokes  polarization  vectors  and  the  4*4  Mueller  matrices 
encountered  in  the  study  of  polarization  in  Chapter  6.   For  a  more  general  treatment  the 
reader  should  consult  mathematics  or  mathematical  physics  texts  [7,18]. 

A  Mueller  matrix  is  a  two-dimensional  array  of  16  numbers,  the  symbol  for  any  one  of 
which  may  be  labeled  by  two  subscripts.   The  first  subscript  denotes  the  row  and  the  second 
the  column  of  the  element  in  question.   We  can  write  the  matrix  f\     in  any  of  the  following 
equivalent  ways: 


A  -  [Arc]   = 


A00  A01  A02  A03 

A10  All  A12  A13 

A  A  A  A 

20  21  22  23 


LA30 


AAA 
31  32  33 


(A6.1) 


The  product   of   two  matrices     /\     and     $      is   a  new  matrix 

C  =  A-B 

each  of  whose  elements  is  computed  by  the  formula: 


(A6.2) 


Crc  =  5>rk-Bkc, 

k=0 


r ,c  =  0,  1,  2,  3. 


(A6.3) 


Unfortunately  our  page  isn't  wide  enough  to  display  the  complete  product  matrix  Q  in  one 
block;  however,  if  we  may  be  permitted  to  break  it  between  the  second  and  third  columns  we 
can  show  it  all. 


C  =  [c  ]  = 

re 


A00'B00  +  A01'B10  +  A02-B20  +  A03-B30 

A10-B00+A11-B10+A12-B20+A13'B30 

A20'B00  +  A2l'B10  +  A22#B20  +  A23#B30 

LA30-B00.  +  A31'B10  +  A32'B20  +  A33'B30 


A00'B01  +  A01'B11  +  A02-B21  +  A03'B31 
A10'B01  +  AirBll  +  A12'B21  +  A13'B31 
A20'B01  +  A2l'Bll  +  A22-B21  +  A23'B31 
A30*B01  +  A3l'Bll  +  A32-B21  +  A33*B31 


(A6.4) 


A00-B02-+  A01'B12  +  A02-B22  +  A03'B32 
A10-B02+All'B12+A12'B22+A13-B32 
A20'B02+A21-B12+A22-B22+A23-B32 
A30*B02  +  A3l'B12  +  A32"B22  +  A33*B32 


A00*B03  +  A01*B13  +  A02*B23  +  A03#B33 
A10'B03  +  All'B13  +  A12-B23  +  A13*B33 
A20'B03  +  A21'B13  +  A22'B23  +  A23'B33 
A30'B03  +  A3l'B13  +  A32*B23  +  A33*B33 
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We  may  think  of  a  Stokes  vector 


L  -   [Lr]  - 


LL3J 

as  a  matrix  in  which  all  columns  after  the  first  column  are  zero.   Viewed  this  way 
product  of  a  Mueller  matrix  by  a  Stokes  vector 


(A6.5) 


the 


L'  =  A'L  (A6.6) 

follows  the  formula  given  above.   Alternatively  we  may  write  a  special  formula  for  this 
product  as 


-'  =  Ea 


k=0 


rk  Lk 


r  =  0,  1,  2,  3 


L0  :   A00'L0  +  A0l'Ll  +  A02'L2  +  A03*L3 

Ll  A10'L0  +  All'Ll  +  A12*L2  +  A13*L3 

L2  A20*L0  +  A2l'Ll  +  A22'L2  +  A23*L3 

L3  "  A30*L0  +  A31'L1  +  A32'L2  +  A33*L3 


(A6.7) 


The  meaning  of   the  product   of  a  matrix  or  vector  by  a   scalar    [which   is  any  non-matrix, 
non-vector   quantity   such   as     \,      t(A),      or      ^(L„  +  L    *cos23  +  L    'Sin2g)]      is   that   each 
element   of    the  matrix  or  vector  must   be  multiplied  by   the   scalar.      For   example, 


•A  =  teA    ] 

n  rc 


00 


JsA 


10 


0  1       hk02       ^A03 


^All 
4l 


hA±2       hAu 


^A20        '  ' 

hA3Q 


-2A22       %A23 


%A31       35A32 


%A   „ 


(A6.8) 


Similarly, 


t(X)-|_  = 


x(X)-L 


0 


tW-Lj 
t(X)-L, 


t(X)-L 


3  _ 


(A6.9) 
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Appendix  7.    The  Measurement  of  Mueller  Transmittance  Matrices 

by  John  B.  Shumaker 

In  this  appendix  we  show  how  Mueller  filter  transmittance  matrices  consisting  of  six- 
teen arbitrary  matrix  elements  may  be  determined  experimentally  using  an  unpolarized  source 
of  radiation  and  a  polarization-indifferent  detector.   The  technique  consists  essentially 
of  making  unpolarized  transmittance  measurements  of  polarization  filters  singly,  in  pairs, 
and  in  triples,  as  a  function  of  the  angular  orientations  of  the  filters.   In  practice,  in 
the  presence  of  experimental  measurement  noise,  difficulties  may  arise  in  the  numerical 
solution  of  the  equations  which  follow,  but  we  ignore  such  problems  here.   The  purpose  of 
this  appendix  is  merely  to  demonstrate  in  principle  the  measurability  of  an  arbitrary 
Mueller  matrix  using  only  an  unpolarized  source  and  a  polarization-indifferent  radiometer 
with  no  previously  characterized  polarizing  optical  components. 

The  procedure  simultaneously  measures  three  polarization  filters.   We  assume  that  each 
is  mounted  in  such  a  way  that  it  can  be  reproducibly  inserted  into  and  removed  from  the 
measurement  radiation  beam  and  that  it  can  be  rotated  about  the  beam  axis  with  its  angular 
orientation  read  from  a  scale  fixed  to  the  filter  mount.   If  the  filters  are  labeled   F,  G, 
and  H  we  take  their  Mueller  transmittance  matrices  to  be  perfectly  general: 


F(o°)  = 


6(o°)  = 


00 


01 


02 


03 


F  F     F     F 

10  11    12    13 

P  "C1           "C1           "C1 

20  21    22    23 

F  F     F     F 

30  31    32    33 


G00  G01  G02  G03 

G10  Gll  G12  G13 

G20  G21  G22  G23 

G30  G31  G32  G33 


(A7.1) 


and  similarly  for  H.   Their  dependence  upon  the  orientation  angles   <f>  ,   <J>  ,   and   <}>,   is 
given  by  eq.  (6.39).   These  angles  are  just  the  orientation  angles  read  from  the  scales 
attached  to  the  filter  mounts.   We  make  no  assumption  about  the  relationship  between  the 
<j)  =  0°  directions  (at  which  the  matrices  written  above  apply)  and  the  polarization  axes  of 
the  filters,  if  any. 

The  simplest  measurements  to  make  are  the  transmittances  of  the  individual  filters  as 
measured  with  our  beam  of  unpolarized  light  and  our  polarization-indifferent  detector.   We 
shall  denote  these  measured  transmittances  by  simply  writing  a  bar  over  the  filter  label: 
F,   G,   and  H.   These  transmittances  give  us  directly 


F  "  Foo'  G  =  Goo>  and  H  =  H00' 


(A7.2) 
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The  next  simplest  set  of  measurements  we  can  make  is  of  the  filters  taken  in  pairs.   Let  us 
consider  the  pair   F  and  H  and  pass  the  measurement  beam  through  in  the  direction  from  F 
to  H.   The  measured  pair  transmittance,  which  we  will  call  HF,   is  just  the   (0,0)   element 
of  the  matrix  product 

H(*h)  •  F(*f). 

That  is 

™  ■  (Hoo-Foo  +  ■WV  +  (Hoi'Fio  +  ^-V-Sf  +  (Hoi'F20  -  "oz'V^hf   (A7'3) 

where  C,  c   =  cos2(d), -d)£)  ,   S,£  =  sin2(d>  -<fc.)  ,   and  the  angular  dependence  comes  from  eq. 
nr         h   r     nr         h   r 

(6.39).   If  this  measurement  is  made  at  a  number  of  angular  differences,   cfv-^*  the  results 

h   f 

can  be  solved  as  simultaneous  equations  or  least-squares  fitted  to  obtain  values  for  the 
coefficients  (appearing  in  parentheses)  of  the  constant  and  cosine  and  sine  factors  in 
eq.  (A7.3).   Let  us  call  these  quantities  HF    HF    and  HF  : 

ffil  =  K00*F00  +  H03'F30  (A7'4a) 

^2  "  HorF10  +  H02'F20  (A7'4b) 

^3  =  H01'F20  "  H02*F10-  (A7'4C) 

We  can  also  reverse  either  filter  and  repeat  the  measurements.   Using  eq.  (6.40)  we  see  that 
we  will  get  two  more  sets  of  equations: 

^1  =  H00'F00  +  H03'F03  (A7'5a) 

52  =    H0l'F01  "  H02'F02  (A7'5b) 

53  =-H0rF02-H02'F01  (A7'5C) 
and 

^1  "  H00'F00  +  H30'F30  (A7'6a) 

^2  =  H10-F10  "  H20-F20  (A7'6b) 

^3  "  H10'F20  +  H20-F10  (A7'6C) 

where  we  have  used  a  prime  (')  to  indicate  a  reversed  filter.   Using  the  other  two  pairs  of 
filters  F,G  and  G,H  we  can  obtain  similar  equations  in  products  of  F  and  G  matrix 
elements  and  products  of   G  and  H  matrix  elements.   Since  the   (0,0)   elements  of   p, 
G,   and   H   are  known  [eq.  (A7.2)],  eq.  (A7.5a)  and  its  analogs  in  FG  and   GH  can  be 
solved  for   F   ,   G  .,   and  H  _.   Likewise  from  eq.  (A7.6a)  and  its  analogs   F   ,   G  Q, 
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and  H  n   can  be  obtained.   The  result  for  filter  G  is 


0  3 


(FG  '  -  F-G)'  (GH  '  -G-H) 
HF'  -  H-F 


(A7.7a) 


3  0 


(F'G   -  F-G)- (G'H   -  G-H) 
HTF1  -  H-F 


(A7.7b) 


The  other  solutions  can  be  obtained  by  cyclic  permutations  of  the  symbols  F,  G,  and  H. 
Notice  that  when  G     (or  G   )   is  calculated  we  may  take  either  the  positive  or  negative 
sign  for  the  square  root.   This  ambiguity  arises  because  we  have  not  defined  in  our  experi- 
ment what  we  mean  by  right-circular  polarization.   We  are  free  to  choose  either  direction 
to  be  associated  with  a  positive  value  for  the  L„   Stokes  component.   If  consistency  with 
an  independent  standard  of  handedness  is  necessary  then  a  filter  or  beam  of  known  handed- 
ness is  required.   Once  one  sign  has  been  chosen  then  the  signs  of  all  the  other   (0,3)   and 
(3,0)   matrix  elements  are  uniquely  defined  by  equations  (A7.4a),  (A7.5a),  (A7.6a)  and  their 
analogs. 

Equations  (A7.5b)  and  (A7.5c)  and  the  similar  equations  written  for  the  other  two 
filter  pairs  constitute  six  equations  in  six  unknowns.   The  solution  of  this  system  for 
filter  G,  for  example,  is 


Yl+  Y; 


0  1 


HF2'2  +  HF  '2 


(A7.8a) 


■,  2   = 
"0  2 


Y1+  Y2 


HF„'2  +  HF  '2 


(A7.8b) 


where 


and 


=  [(GH'2  +  GH  '2)-(HF'2  +  HF  '2)'(FG„'2  +  FG _'2)]^ 


Y2  =  HF3'-(GH2' -FG3'  +  GH^  -FG2')+  HF2' •  (GH^ -FG^  -  GH^-FGp 


The  solutions  for   F   and  H  can  be  obtained  by  permuting  the  symbols.   The  solution  of 
the  six  equations  like  eqs.  (A7.6b)  and  (A7.6c)  for  the   (1,0)   and   (2,0)   matrix  elements 
follows  immediately  by  symmetry.   Again  we  see  that  the  sign  of   G   ,  for  example,  can  be 
chosen  arbitrarily  after  taking  the  square  root  of  eq.  (A7.8a).   As  with  G    we  are  free 
to  take  either  sign  but  not  if  we  want  to  be  consistent  with  other  polarization  optics.   If, 
for  example,   G   is  primarily  a  linear  polarizer  which  is  more  transparent  to  horizontally 
polarized  light  when   4>   =  0   than  to  vertically  polarized  light  then  G    should  be  chosen 
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positive.   A  suitable  source  of  partially  horizontally  polarized  light  is  light  specularly 
reflected  from  a  glass  surface  oriented  so  as  to  produce  a  vertical  plane  of  incidence. 
Similarly,  if  at   <j>   =0   it  is  more  transparent  to  +45°  polarized  light  than  to  +135° 
polarized  light  then  G    should  be  taken  positive.   Once  these  two  signs  are  established 
all  the  others  are  determined  by  eqs.  (A7.4),  (A7.5),  (A7.6),  etc. 

Finally  we  measure  the  transmittance  of  all  three  filters  together.   Let  us  assume 
that  the  beam  passes  through  them  in  the  order   F,G,H   so  that  the  measured  transmittance, 
HGF   is  given  by  the   (0,0)   element  of  the  matrix  product 

HUh)-G(*g)-F(*f)r 

If  we  perform  the  matrix  multiplication  we  find  that  this   (0,0)   element  is  given  by 


HGF  =  HGF1  +  HGF  -C  ,  +  HGF  -S  ,  +  HGF  '(I   +  HGF  -S, 
1      2   gf      3   gf      4   hg      5   hg 

+  HGF.-C   -C,   +  HGF  -C  . • Su   +  HGF-S   -C.   +  HGF  -S   -S, 
6  gf  hg      7   gf  hg      8  gf  hg      9  gf  h£ 


(A7.9) 


where 


HGF1  =  H00-G00'F00  +  H00'G03-F30  +  H03'G30*F00  +  H03-G33'F30 

^2  "  H00*G01-F10  +  H00-G02'F20  +  H03'G31-F10  +  H03'G32-F20 

^3  "  H00'G0l'F20  "  H00'G02-F10  +  "oS^Sl'^O  "  H03'G32'F10 
HGF4  =   H01-G10-F0()  +  H01-G13-F30  +  H^-G^-F^  +  H^-G^-F^ 

^5   =  -H02-G10-F00  "  H02-G13'F30  +  ^/So^OO  +  H0rG23'F30  (A7'10) 

HGF6  =   H01-Gu-F10  +  HQ1-G12-F20  +  H^-G^-F^  +  H^-G^-F^ 

HGF7  =  "H02*Gll"F10  "  H02'G12'F20  +  H01*G2l'F10  +  H0l"G22'F20 

HGF8  =  H0l'Gll'F20  "  H0l'G12'F10  +  H02*G2l'F20  "  H02*G22'F10 
HGF9  =  -H02-Gn-F20  +  HQ2-G12-F10  +  H^-G^-F^  -  HQ1-G22-F10 

and   C   =  cos2((()  -§   ) ,   S   =  sin2(<J)  -$   ) ,   etc.   Again  we  assume  that  by  least  squares 
fitting  or  by  solution  of  simultaneous  equations  we  have  obtained  numerical  values  for  the 
coefficients   HGF  ,   HGF„,   etc.,  of  all  the  trigonometric  factors  in  eq.  (A7.9);  that  is, 
we  assume  that  we  know  the  left-hand  sides  of  eqs.  (A7.10).   In  the  first  of  eqs.  (A7.10) 
we  now  know  everything  except   G    so  it  can  be  solved  for  directly.   The  second  and  third 
of  eqs.  (A7.10)  are  a  pair  of  equations  in  which  the  only  remaining  unknowns  are  G    and 
G„_.   Likewise  the  next  two  equations  contain  only  two  unknowns,   G    and   G   .   The  last 
four  equations  of  the  set  (A7.10)  contain  the  four  unknowns   G   ,  G  „,  G   ,  and  G   .   Thus, 
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all  the  rest  of  the  matrix  elements  of  the  Mueller  matrix  for  the  filter  G  can  be  deter- 
mined from  these  nine  coefficient  values.   Measurements  of  the  three-filter  train  in  the 
orders  G,  H,  F  and   H,  F,  G  will  permit  the  remaining  matrix  elements  of  the  other  two 
filters  to  be  determined.   We  won't  attempt  to  write  out  explicit  solutions  for  the  unknowns 
appearing  in  the  set  of  equations  (A7.10).   The  procedure  to  obtain  the  solutions  is 
obvious  and  straightforward,  but  the  entire  data  reduction  process  for  this  general  case  is 
clearly  best  left  to  a  computer,  if  possible. 

We  have  assumed  that  the  matrix  elements  in  the  first  row  and  first  column  are  non-zero. 
If  zeros  appear  here  some  of  the  matrix  elements  may  become  indeterminate.   For  example,  in 
eq.  (A7.7a)  if  either  F    or  H    vanishes  the  denominator  and  one  factor  of  the  numera- 
tor will  both  vanish  and   G    cannot  be  evaluated.   The  physical  reason  for  this  is  clear. 
The  experimental  procedure,  in  effect,  requires  generating  partially  circularly  polarized 
light  with  one  filter  and  then  analyzing  it  with  a  second  filter.   If  the  first  filter  fails 
to  produce  circularly  polarized  light,  no  information  can  be  obtained  about  the  circular 
polarization  behavior  of  the  second  filter.   Similar  arguments  apply  to  the  other  first 
column  and  first  row  matrix  elements.   The  extreme  case  occurs  when  all  the  first  column 
and  first  row  elements  of  one  filter,  say  G,   vanish.   Then  G   is  totally  opaque  and 
obviously  contributes  nothing  to  the  measurement  of  filters   F  and   H.   With  only  two 
filters,   F  _   and  H  „   can  be  measured  as  before,  but  equations  (A7.4),  (A7.5),  and  A7.6) 
provide  only  9  equations  in  12  unknowns  and  contain  no  information  about  the  18  other  matrix 
elements  which  do  not  lie  in  the  first  row  or  column.   When  the  equations  become  insoluble 
because  of  the  occurrence  of  zero  denominators  one  must  either  abandon  the  attempt  to 
measure  the  inaccessible  matrix  elements  or  find  other  filters  which  don't  have  such  limi- 
tations. 

When  the  matrix  elements  for  an  optical  component  have  been  determined  as  outlined 
above  they  can  be  inserted  in  eq.  (6.39)  to  obtain  the  orientation  dependence  of  the  Mueller 
matrix.   In  general,  even  for  ideal  optical  components  the  result  will  not  agree  with,  for 
example  eq.  (6.15),  but  this  is  only  because  the  angle  between  the  polarization  axis  of  the 
device  and  the   <j>  =  0°  direction  was  assumed  arbitrary.   The  direction  of  the  polarization 
axis  can  be  found,  if  desired,  by  finding  the  angle   cf>°  at  which  the  experimental  Mueller 
matrix  most  nearly  resembles  the  matrix  for  the  ideal  component.   For  example,  for  an  ideal 
linear  polarizer  at   <J>  =  0°  the   (0,2)   element  vanishes  so  for  a  real  polarizer  we  must 
find  the  angle   <j>°  at  which  the   (0,2)   element  [from  eq.  (6.39)]  of  the  experimental 
matrix  vanishes.   That  is, 

G      •sin2<))0   +  G      -cos2cf>0    =   0 

or 

<f>°   =   -Jsarctan(G02/G01). 

This  is  the  direction  of  the  polarization  axis  and  if  the  orientation  angle  for  an  ideal 
component  is  measured  from  this  direction  instead  of  from  <f>  =  0°  its  experimentally 
determined  Mueller  matrix  should  agree  with  that  given  by  eq.  (6.15). 
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organizations,  and  regulatory  bodies. 
Special  Publications — Include  proceedings  of  conferences 
sponsored  by  NBS,  NBS  annual  reports,  and  other 
special  publications  appropriate  to  this  grouping  such 
as  wall  charts,  pocket  cards,  and  bibliographies. 
Applied  Mathematics  Series — Mathematical  tables,  man- 
uals, and  studies  of  special  interest  to  physicists,  engi- 
neers, chemists,  biologists,  mathematicians,  com- 
puter programmers,  and  others  engaged  in  scientific 
and  technical  work. 

National  Standard  Reference  Data  Series — Provides 
quantitative  data  on  the  physical  and  chemical  proper- 
ties of  materials,  compiled  from  the  world's  literature 
and  critically  evaluated.  Developed  under  a  world-wide 
program  coordinated  by  NBS.  Program  under  authority 
of  National  Standard  Data  Act  (Public  Law  90-396). 


NOTE:  At  present  the  principal  publication  outlet  for 
these  data  is  the  Journal  of  Physical  and  Chemical 
Reference  Data  (JPCRD)  published  quarterly  for  NBS 
by  the  American  Chemical  Society  (ACS)  and  the  Amer- 
ican Institute  of  Physics  (AIP).  Subscriptions,  reprints, 
and  supplements  available  from  ACS,  1155  Sixteenth 
St.  N.W.,  Wash.  D.  C.  20056. 

Building  Science  Series — Disseminates  technical  infor- 
mation developed  at  the  Bureau  on  building  materials, 
components,  systems,  and  whole  structures.  The  series 
presents  research  results,  test  methods,  and  perform- 
ance criteria  related  to  the  structural  and  environmental 
functions  and  the  durability  and  safety  characteristics 
of  building  elements  and  systems. 

Technical  Notes — Studies  or  reports  which  are  complete 
in  themselves  but  restrictive  in  their  treatment  of  a 
subject.  Analogous  to  monographs  but  not  so  compre- 
hensive in  scope  or  definitive  in  treatment  of  the  sub- 
ject area.  Often  serve  as  a  vehicle  for  final  reports  of 
work  performed  at  NBS  under  the  sponsorship  of  other 
government  agencies. 

Voluntary  Product  Standards — Developed  under  proce- 
dures published  by  the  Department  of  Commerce  in  Part 
10,  Title  15,  of  the  Code  of  Federal  Regulations.  The 
purpose  of  the  standards  is  to  establish  nationally  rec- 
ognized requirements  for  products,  and  to  provide  all 
concerned  interests  with  a  basis  for  common  under- 
standing of  the  characteristics  of  the  products.  NBS 
administers  this  program  as  a  supplement  to  the  activi- 
ties of  the  private  sector  standardizing  organizations. 
Consumer  Information  Series — Practical  information, 
based  on  NBS  research  and  experience,  covering  areas 
of  interest  to  the  consumer.  Easily  understandable  lang- 
uage and  illustrations  provide  useful  background  knowl- 
edge for  shopping  in  today's  technological  marketplace. 

Order  above  NBS  publications  from:  Superintendent 
of  Documents,  Government  Printing  Office,  Washington, 
D.C.  201,02. 

Order  following  NBS  publications— NBSIR's  and  FIPS 
from  the  National  Technical  Information  Services, 
Springfield,  Va.  22161. 

Federal  Information  Processing  Standards  Publications 
(FIPS  PUBS) — Publications  in  this  series  collectively 
constitute  the  Federal  Information  Processing  Stand- 
ards Register.  Register  serves  as  the  official  source  of 
information  in  the  Federal  Government  regarding  stand- 
ards issued  by  NBS  pursuant  to  the  Federal  Property 
and  Administrative  Services  Act  of  1949  as  amended, 
Public  Law  89-306  (79  Stat.  1127),  and  as  implemented 
by  Executive  Order  11717  (38  FR  12315,  dated  May  11, 
1973)  and  Part  6  of  Title  15  CFR  (Code  of  Federal 
Regulations). 

NBS  Interagency  Reports  (NBSIR) — A  special  series  of 
interim  or  final  reports  on  work  performed  by  NBS  for 
outside  sponsors  (both  government  and  non-govern- 
ment). In  general,  initial  distribution  is  handled  by  the 
sponsor;  public  distribution  is  by  the  National  Techni- 
cal Information  Services  (Springfield,  Va.  22161)  in 
paper  copy  or  microfiche  form. 


BIBLIOGRAPHIC  SUBSCRIPTION  SERVICES 


The  following  current-awareness  and  literature-survey 
bibliographies  are  issued  periodically  by  the  Bureau: 
Cryogenic  Data  Center  Current  Awareness  Service.  A 

literature   survey  issued  biweekly.  Annual   subscrip- 
tion: Domestic,  $25.00  ;  Foreign,  830.00'. 
Liquified  Natural  Gas.  A  literature  survey  issued  quar- 
terly. Annual  subscription:   $20.00. 


Superconducting  Devices  and  Materials.  A  literature 
survey  issued  quarterly.  Annual  subscription:  $30.00  . 
Send  subscription  orders  and  remittances  for  the  pre- 
ceding bibliographic  services  to  National  Bureau  of 
Standards,  Cryogenic  Data  Center  (275.02)  Boulder, 
Colorado  80302. 


U.S.  DEPARTMENT  OF  COMMERCE 
National  Bureau  of  Standards 

Washington.  DC-  20234 


OFFICIAL  BUSINESS 

Penalty  for  Private  Use.  S300 


POSTAGE  AND  FEES  PAIO 

U.S.  DEPARTMENT  OF  COMMERCE 

COM-215 


SPECIAL  FOURTH-CLASS  RATE 
BOOK 


